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A  CLOUD  CHAMBER  STUDY  OF  THE  EFFECT  WHICH  NONPRECIPITATING 
CLOUDS  HAVE  ON  THE  AEROSOL  SIZE  DISTRIBUTION 


1.  INTRODUCTION 

It  is  well  known  that  background  aerosol  affects  cloud  microstructure,  and  hence  influences  cloud 
optical  properties  and  radiative  transfer,  a  subject  of  great  interest  in  climate  research.  Cloud  droplets 
form  on  preexisting  aerosol  particles  called  cloud  condensation  nuclei  (CCN).  Whether  or  not  an  aerosol 
particle  will  act  as  a  CCN  in  a  given  cloud  depends  in  a  complicated  way  on  the  cooling  rate,  and  the 
number,  size  distribution,  and  composition  of  the  particles.  For  the  composition  of  particles  typically 
found  in  the  atmosphere,  particles  larger  than  about  0.03  to  0.08  jtm  radius  (depending  on  the  maximum 
cloud  supersaturation)  will  be  activated  in  typical  clouds. 

Less  understood  and  less  recognized  is  the  fact  that  nonprecipitating  cloud  cycles  are  a  significant 
factor  in  shaping  the  background  aerosol  size  distribution.  More  than  90%  of  the  clouds  in  the 
atmosphere  (Pruppacher  1986)  reevaporate  without  the  formation  of  precipitation.  This  means  that  a 
parcel  of  air  will  be  processed  through  many  nonprecipitating  clouds  (on  average)  before  it  is 
incorporated  into  a  precipitating  cloud  where  the  aerosol  can  be  removed  by  precipitation  scavenging.  The 
cloud  droplet  radius  (mass)  is  approximately  100  (106)  times  larger  than  that  of  the  "dry"  CCN.  During 
the  cloud  phase,  the  droplet  can  absorb  soluble  trace  gases  (e.g.,  sulfur  dioxide  and  oxides  of  nitrogen) 
and  oxidants  (e.g.,  ozone  and  hydrogen  peroxide)  which  react  in  the  liquid  phase  to  form  less  volatile 
products  that  remain  as  residue  after  the  cloud  evaporates,  thus  causing  an  increase  in  aerosol  mass.  In 
an  early  paper,  Hobbs  (1971)  suggested  that  the  production  of  sulfates  in  cloud  droplets,  followed  by  the 
release  of  these  particles  when  the  droplet  evaporates,  could  be  an  important  mechanism  for  increasing 
the  efficiency  of  CCN.  In  a  subsequent  paper,  Hegg  et  al.  (1980)  found  an  increase  in  the  concentrations 
of  CCN  downwind  of  a  wave  cloud  and  an  increase  in  the  sulfate  load  of  about  0.5  ng  m‘3  from  filter 
samples.  Measurements  of  the  particle  size  distribution  at  radii  greater  than  about  0.25  /im  with  an  optical 
particle  counter  suggested  that  particles  of  about  0.5  ftm  were  responsible  for  the  increase.  However, 
the  size  distribution  data  was  much  less  convincing  than  the  CCN  and  filter  loading  data  and  probably 
did  not  go  to  sufficiently  small  sizes  to  see  the  region  where  most  conversion  occurs.  In  addition  to 
conversion  within  the  cloud  droplets,  smaller  unactivated  particles  (interstitial  aerosol)  will  diffuse  to  the 
large  surface  area  provided  by  the  cloud  droplets.  The  mass  of  the  smaller  particle  is  incorporated  into 
the  CCN  after  evaporation.  The  former  process  (i.e.,  a  liquid  phase  reaction)  has  a  much  larger  effect 
on  the  mass  of  CCN  than  does  the  latter. 

The  growth  of  CCN  mass  during  nonprecipitating  cloud  cycles  has  a  tendency  to  produce  a 
separation  in  the  size  distribution  at  the  critical  particle  radius  which  corresponds  to  the  supersaturation 
in  the  cloud.  Hoppel  et  al.  (1983,  1986,  1990)  hypothesized  that  this  mechanism  caused  the 
double-peaked  size  distribution  which  they  observed  in  the  0.01  to  0.2  j*m  size  range  in  the  marine 
boundary  layer.  The  minimum  between  the  two  peaks  is  between  0.06  to  0. 1  jim  radius.  This  persistent 
double-peaked  feature  is  observed  both  over  the  Atlantic  (Hoppel  et  al.  1990)  and  Pacific  Oceans  (Hoppel 
and  Frick  1990).  While  the  double-peaked  feature  is  not  as  conspicuous  over  land,  we  believe 
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nonprecipitating  clouds  play  an  equal,  if  not  more  important,  role  in  supplying  aerosol  mass  over 
continental  regions.  The  minimum  between  the  two  peaks  is  obscured  by  locally  distributed  sources  over 
continental  regions  and  the  fact  that  the  relationship  between  size  and  critical  supersaturation  is  less  well 
defined  for  continental  aerosol  whose  composition  is  less  well  defined.  However,  there  is  often  a  large 
amount  of  mass  created  in  the  0.1  to  0.5  /*m  radius  range  over  continents  which,  under  some 
circumstances,  is  best  explained  in  terms  of  the  nonprecipitating  cloud  process. 

This  report  describes  experiments  carried  out  in  Calspan’s  590  m3  chamber  to  study  the  effect  of 
nonprecipitating  clouds  on  the  aerosol  size  distribution  under  well -controlled  conditions.  The  use  of  a  very 
large  chamber  is  required  because  cloud  droplets  which  are  of  the  order  of  10  n m  will  fall  about  3  m  in 
5  minutes.  The  change  in  the  aerosol  size  distribution  during  a  cloud  cycle  which  resulted  from  the  liquid 
phase  oxidation  of  S02  by  ozone  in  cloud  droplets  was  measured  and  the  results  compared  to  those 
obtained  by  modeling.  Oxidation  of  SC^  by  hydrogen  peroxide  during  a  cloud  cycle  was  also  observed 
but  a  detailed  quantitative  comparison  with  theory  is  not  possible  because  no  capability  was  available  to 
measure  hydrogen  peroxide  concentrations. 

2.  AEROSOL  GROWTH  LAWS  FOR  LIQUID  PHASE  CONVERSION  OF  SULFUR  DIOXIDE 

TO  SULFATE  IN  CLOUD  DROPLETS 

From  the  rate  at  which  trace  gases  and  oxidants  react  in  the  cloud  droplet  to  form  involatile  products, 
growth  laws  for  the  equivalent  "dry"  size  can  be  derived.  This  can  be  done  provided  the  following  are 
known: 


•  the  liquid  phase  reaction  rates, 

•  the  gas  phase  concentrations  of  the  reactants, 

•  the  Henry’s  law  constants, 

•  the  cloud  droplet  size,  and 

•  the  composition  of  the  CCN. 

2.1  Growth  Law  for  the  Oxidation  of  Sulfur  Dioxide  by  Ozone  in  a  Cloud  Droplets 

The  rate  of  oxidation  of  S02  by  03  in  bulk  water  which  is  in  Henry’s  law  equilibrium  with  SO^  and 
03  in  the  surrounding  air  is  given  by  Hoffmann  (1986)  as 


+  k 


KxK2 

[*f 


H0HS  P0rs 


0) 


where  [S04]  is  the  concentration  of  sulfate  in  moles  l'1,  [H+]  is  the  hydrogen  ion  concentratirn  in  moles 
I1,  H0  and  Hs  are  the  Henry’s  law  constants  (moles  1 1  atm'1);  P0  and  Ps  are  the  ambient  vapor  pressures 
for  ozone  and  sulfur  dioxide  in  atmospheres;  k$,  kx,  and  are  the  rate  constants  for  reaction  of  03  with 
S02  •  H20,  HS03‘,  and  S032",  respectively;  and  Kx  and  K2  are  the  first  and  second  dissociation  constants 
for  dissolved  S02  in  moles  l'1. 


In  terms  of  the  total  conversion  rate  in  a  droplet  of  radius  R  (in  cm),  Eq.  (1)  becomes 
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where  [S04]  is  the  moles  of  sulfate  in  the  droplet.  In  Eq.  (2),  values  of  the  constants  for  T  =  25°  C 
given  by  Hoffmann  (1986)  and  Jacob  and  Hoffmann  (1983)  have  been  adopted  and  the  first  term  on  the 
right  of  Eq.  (1)  has  been  dropped  because  it  has  a  negligible  effect  on  the  growth  rates.  Assuming 
complete  ionization  of  sulfuric  acid  (H2S04  -*  H+  +  HS04),  the  hydrogen  ion  concentration  and  the 
moles  of  sulfate  in  the  droplet  can  be  expressed  in  terms  of  the  mass  of  sulfuric  acid  Ma  as 


and 


[H*]  *  2.436  — 
1  J  R3 


(3) 
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where  Mw  is  the  molecular  weight  of  sulfuric  acid.  Equation  (3)  is  valid  only  when  the  pH  is  less  than 
about  6.5.  The  mass  growth  rate  of  sulfuric  acid  in  a  cloud  droplet  can  then  be  expressed  as 
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Equation  (5)  can  be  integrated  so  that  the  acidic  mass  at  a  later  time  can  be  found  for  a  droplet  of 
fixed  radius  in  an  atmosphere  of  known  S02  and  03  concentrations. 


The  two  terms  on  the  right  of  Eqs.  (2)  and  (5)  are  equal  when  the  pH  of  the  droplet  is  about  3.6. 
For  most  calculations,  only  the  [H+]'2  term  is  important.  This  is  because  at  low  pH  (high  [//+] 
concentrations),  the  growth  is  negligible  on  the  time  scale  of  interest  in  this  study.  For  higher  pH  where 
growth  rates  are  significant,  the  second  term  makes  the  dominant  contribution.  If  only  the  second  term 
is  retained,  integration  yields 


M«(t)  =  M]( 0)  +  3.92  x  1( r3R9P0Pst  <6) 

where  Ma( 0)  is  the  initial  mass  of  the  sulfuric  acid  in  the  droplet. 

Since  sulfuric  acid  is  hygroscopic,  even  at  low  relative  humidities  (RH),  it  is  convenient  to  write  the 
initial  mass  of  sulfuric  acid  in  terms  of  the  total  initial  mass  MT(Q)  (at  some  reference  relative  humidity) 
as 


Ma(0)  =  ceMT(  0)  (7) 

where  c  is  the  weight  fraction  of  the  sulfuric  acid  in  solution  at  the  reference  RH,  and  e  the  fraction  of 
acidic  solution  in  the  particle.  Thus,  (l-e)MT(0)  represents  the  fraction  of  nonacidic  material  in  the 
original  CCN. 

Likewise,  if  MT(t)  is  the  total  "dry"  mass  of  the  particle  at  time  t  (i.e.,  the  mass  remaining  in  the 
droplet  after  evaporating  to  its  equilibrium  size  at  the  reference  humidity),  then 

Ma(t)  =  c[MT(t)  -  (1  -f)Mr(0)]  (8) 
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where,  for  lack  of  a  better  assumption,  we  have  assumed  that  the  nonacidic  fraction  has  the  same 
hygroscopicity  as  does  the  sulfuric  acid. 

The  reference  relative  humidity  is  taken  to  be  the  relative  humidity  in  the  differential  mobility  size 
spectrometer.  Equation  (6)  can  then  be  written  as 


[cMt(T)  -  c(l-e)Mr(0)]3  =  (ceAf^O)]3  +  3.92  x  10 ^R9P0Pst.  (9) 

If  the  original  CCN  were  totally  acidic  (e  =  1),  Eq.  (9)  could  be  written  in  terms  of  the  radius  as 
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where  r  is  the  "dry"  radius  and  p4  is  the  density  of  the  acid  at  the  reference  RH.  The  radial  growth  law 
can  then  be  written  as 


dr  =  1,31  x  IQ'3  R » 

dt  3  (4/3  v  pa  cf  r8 


(11) 


It  is  evident  that  the  growth  rate  increases  very  rapidly  with  increasing  cloud  droplet  radius  (R9) 

and  decreases  rapidly  with  increasing  particle  radius  (r~8).  This  decrease  with  increasing  particle  radius 
is  much  more  rapid  than  any  of  the  proposed  gas-phase  growth  laws  (McMurry  and  Wilson  1982).  (The 

most  rapidly  decreasing  gas-phase  growth  law  is  for  difftisional  growth  and  is  proportional  to  r'1 .)  Once 
the  growth  rate  reaches  a  threshold  where  the  rate  is  observable  on  the  time  scale  of  the  experiment, 
further  increase  of  R  dramatically  increases  the  growth  rate  of  r,  so  that  r  then  tracks  the  increase  in  the 
cloud  droplet  radius  R.  Physically,  this  tracking  is  largely  caused  by  the  pH  effect  of  the  reaction.  As 
the  cloud  droplet  radius  increases,  dilution  occurs,  increasing  the  pH;  the  reaction  rate  increases  again, 
lowering  the  pH  and  slowing  the  conversion  rate. 

This  tracking  is  better  illustrated  by  Eq.(10).  If  there  is  any  significant  growth,  then  r9(t)>  >r( 0)9 
and  Eq.  (10)  becomes 


KO  _ 
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Values  of  c  =  0.5  and  pt  =  1.4,  corresponding  to  our  reference  humidity  of  34%  (humidity  in  the 
DMA  size  spectrometer),  have  been  used  in  Eq.  (12). 


The  dependence  of  the  ratio  r(t)/R  on  P O’  Ps,  and  time  in  the  cloud  is  very  weak.  For  conditions 
typical  of  those  in  our  experiments  ( P0  =  90  ppb,  Ps  =  6  ppb,  and  t  =  240  s),  this  ratio  becomes 


—  =  0.137.  (13) 
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Changing  t  or  their  product  by  an  order  of  magnitude  changes  this  ratio  by  less  than  30% . 

This  ratio  corresponds  to  a  pH  in  the  droplet  of  about  4.5.  For  5  to  10  p. m  cloud  droplets,  this 
corresponds  to  dry  radii  of  0.07  to  0. 14  p m. 

Figure  1  shows  solutions  to  Eq.  (10)  for  three  cloud  droplet  radii  ( R  =  4,  6,  and  8  pm).  The  time 
in  cloud  is  4  min.  Nuclei  larger  than  about  0. 15  pm  do  not  grow  even  in  the  8  /4m  radius  cloud  droplet 
because  the  low  pH  caused  by  the  original  nucleus  is  sufficient  to  turn  off  the  conversion  of  S02.  Drops 
which  form  on  small  nuclei  are  sufficiently  dilute  (high  pH)  that  conversion  of  S02  can  occur.  Figure 
2  shows  the  growth  of  the  "dry"  size  as  a  function  of  time  as  calculated  from  Eq.  (10),  starting  with  a 
nucleus  whose  radius  is  0.04  pm.  The  growth  is  initially  very  rapid  in  the  larger  droplet  but  is  quickly 
diminished  as  the  pH  of  the  droplet  is  decreased  by  conversion  of  SOj  to  sulfuric  acid. 

The  above  arguments  are  for  acidic  nuclei.  For  nonacidic  nuclei,  acidic  mass  can  be  formed  before 
the  pH  is  sufficiently  low  to  quench  the  conversion.  For  the  case  of  a  nonacidic  nucleus  (e  =  0),  Eq. 
(9)  becomes 


r(t) 

3  -  fr(0)l 

3 

+  • 

3  92  * 10'3  p0ps  , 

R 

l  R 

(4/3  7T  Pacf 

(14) 


where  the  density  p0  of  the  nonacidic  material  is  assumed  to  be  the  same  as  that  of  the  sulfuric  acid  at 
the  reference  relative  humidity  and  to  have  the  same  concentration  of  associated  water.  Figure  3  shows 
the  change  in  the  "dry"  radius  before  and  after  a  cloud  for  the  same  conditions  as  in  Fig.  1  except  that 
here  the  starting  nuclei  are  not  acidic. 

2.2  Growth  Law  for  Oxidation  of  Sulfur  Dioxide  by  Hydrogen  Peroxide  in  a  Cloud  Droplet 

The  rate  of  sulfate  formation  due  to  aqueous-phase  oxidation  of  SO2  by  H202  can  be  written 
(Overton  1985)  as 


-  k  [ff2o2]  [s(/v>) 


(15) 


where  S(IV )  is  the  sulfur  dioxide  in  the  water  as  hydrated  sulfur  dioxide,  bisulfite  ion,  and  sulfite  ion; 
and  k  is  the  effective  second-order  rate  constant  for  the  H202  -  S(IV)  reaction.  Since 


we  have 


[S(IV)\  =  [S02(aq)} 


2 

[H  f 


dt 


K  K,K , 

+  _ + 

[/r]  [h+]2 


s^h 


(16) 
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Fig.  1  —  Growth  of  "equivalent  dry  radius"  of  nucleus  during  a  4 
min  cloud  cycle  for  cloud  droplets  of  4,  6,  and  8  /tm  radius.  S02 
and  03  concentrations  are  6  ppb  and  90  ppb,  respectively.  The 
equivalent  dry  radius  is  taken  to  be  the  radius  of  a  sulfuric  acid 
nucleus  at  34%  RH. 


Fig.  2  —  Growth  of  equivalent  dry  radius  as  a  function  of  time  in  cloud  droplets  of 
4,  6,  and  8  >im  radius.  Conditions  are  the  same  as  in  Fig.  1. 
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Fig.  3  —  Same  as  Fig.  1  except  here  the  starting  nucleus 
is  assumed  to  be  nonacidic 


where  Hh  and  Hs  are  Henry’s  Law  constants  for  H202  and  S02,  respectively,  and  Ph  and  Ps  are  the 
environmental  pressures  of  H202  and  S02  in  atmospheres.  If  we  adapt  Overton’s  (1985)  expression  for 
k,  the  product  of  k  and  (1  +Kl/[H+]+KlK2/[H+]2)  has  a  nearly  constant  value  of  7.7  x  105  s'1  atm'1  for 
pH  values  greater  than  about  2.0.  Evaluating  Eq.  (17)  and  using  Hh  =  9.7  x  104  moles  l'1  atm'1 
(Chameides  1984)  and  Hs  =  1.24  moles  l'1  atm'1  gives 

=  9.3  x  1010  PHPS.  (lg) 

dt  H  s 

Converting  Eq.  (18)  from  moles  per  liter  to  moles  per  droplet  where  the  droplet  has  a  radius  of  R 
gives 


d[S04]  _  4/3  tt 
dt  1000 


(9.3  x  1010)  R3  PHPS. 


(19) 


Using  Eq.  (4)  to  convert  moles  of  sulfate  to  total  mass  of  sulfuric  acid  Ma  in  the  droplet,  Eq.  (19) 
becomes 


dMa 

dt 


=  3.9  x  1010  R 3  PHPS. 


(20) 
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Since  the  conversion  rate  is  proportional  to  R3,  the  conversion  of  S02  by  H2C>2  is  expected  to  be 
proportional  to  the  liquid  water  content  of  the  cloud  provided  the  reactants  are  not  depleted;  this  was  not 
the  case  for  liquid-phase  conversion  of  S02  by  ozone.  Integration  of  Eq.  (20)  yields 

Ma(t)  =  Ma(o)  +  3.9  x  1010  R3  PHPS  t.  (21) 


Converting  acidic  mass  in  Eq.  (21)  to  total  "dry"  mass  by  using  Eq.(8),  we  obtain 

MT{t)  =  MT(o)  *  -9  *  l?10  R 3  PHPS  t. 


(22) 


In  terms  of  radius 


r(t)  . 
R 


+ 


3.9  x  IQ10 
4/3  7T  pac 


1/3 


(23) 


where,  for  simplicity,  we  have  assumed  that  the  nonacidic  component  in  the  original  nuclei  has 
approximately  the  same  density  as  the  acid  at  the  reference  humidity.  In  this  case,  the  radial  growth  law 
is  given  by 

dr  ,  3.9  X  I0‘°  £  po  (24) 

dt  4  tt  pac  r 2 


The  growth  rate  is  proportional  to  the  volume  of  the  cloud  droplet  and  inversely  proportional  to  r2. 

Figure  4  shows  the  growth  of  the  equivalent  dry  radius  at  the  reference  RH  over  a  period  of  4  min 
in  cloud  droplets  of  4,  6,  and  8  pm  radius,  similar  to  that  shown  in  Fig.  2  for  03  oxidation.  For  the 
H202  oxidation,  the  same  SOj  concentration  of  6  ppb  was  used  but  the  H202  concentration  was  assumed 
to  be  only  1  ppb.  Still,  the  growth  for  the  H202  oxidation  of  S02  is  much  faster  than  with  90  ppb  of 
ozone. 


TIME  (min) 

Fig.  4  —  Growth  of  equivalent  dry  radius  in  cloud  droplets  of  4,  6,  and  8  fitn  radius  for 
the  case  of  6  ppb  S02  and  1  ppb  H202 
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Figure  5  compares  the  radius  of  the  CCN  before  and  after  a  4  min  cloud  cycle  for  the  conditions 
stated  above  ( PH  =  lppb).  Here  the  growth  is  not  limited  by  the  pH. 


RADIUS  BEFORE  CLOUD  (pm) 

Fig.  5  —  Growth  of  equivalent  dry  radius  of  nucleus  during  a  4  min 
cloud  cycle  for  oxidation  by  H202.  S02  and  H202  concentrations 
were  taken  to  be  6  ppb  and  1  ppb,  respectively. 


If  r*  >  >  r03,  then  the  equivalent  dry  radius  is  proportional  to  the  cloud  drop  radius  and  to  the  cube 
root  of  the  time  (Eq.  (23)).  From  Eqs.  (3)  and  (21)  it  is  obvious  that  the  pH  in  the  droplet  is  inversely 
proportional  to  the  time  and  independent  of  the  radius  of  the  cloud  droplet;  hence, 

[/T]  =  9.5  x  10'°  PHPst. 

For  the  conditions  cited  and  4  min  in  a  cloud  the  pH  is  about  3.9.  (When  pH  drops  below  about 
3.0,  our  original  assumption  that  the  product  k(l  +  K{/[H+]+ KlKll[H+]2)  is  constant  would  not  be  valid.) 

It  must  be  remembered  that  the  assumed  conditions  of  a  cloud  droplet  growing  in  an  environment 
of  constant  reactant  gases  is  unrealistic  because  the  1  ppb  of  H202  would  rapidly  be  used  up  and  the 
conversion  would  cease  (by  the  time  about  4  pg  m"3  of  sulfate  aerosol  were  formed) 

2.3  Depletion  of  Reactants 


The  above  analysis  assumes  constant  concentrations  of  reactants.  During  the  conversion  process,  the 
reactants  are  depleted;  for  example,  the  conversion  of  1  ppb  of  SO^  will  produce  about  4.4  fig  m 3  of 
sulfuric  acid.  For  the  case  of  S02  -  03,  the  rate  of  growth  depends  only  weakly  on  the  concentrations 
of  the  reactants  (see  Eq.  (11))  until  one  of  the  reactants  is  nearly  depleted.  If  one  mole  of  gaseous 
reactant  j  produces  one  mole  of  sulf;  2  of  molecular  weight  Mw,  then  the  decrease  in  vapor  pressure  Pj 
(in  atmospheres)  will  produce  an  increase  in  sulfate  mass  density  p,  (in  g  cm"3)  given  by 


dPj- 


2.24  x  104 
M... 


</pr 


(25) 
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For  example,  1  ppb  change  in  SOj  would  produce  approximately  4.4  ng  m'3  of  sulfuric  acid.  In 
terms  of  the  mass  of  a  single  particle  Ma 

dPj  _  _  2.24  x  104  7  dMa  nf^ 

dt  Mw  dt 

where  Z  is  the  number  density  of  cloud  droplets.  If  depletion  is  significant,  then  Eq.  (26)  must  be 
coupled  with  Eq.  (5)  for  S02-03  conversion  or  with  Eq.  (17)  for  S02-H202  conversion.  Depletion  of 
reactants  due  to  in-cloud  conversion  in  our  chamber  experiments  were  not  important  for  the  cases  of  SQ 
oxidation  of  ozone.  However,  in  the  real  atmosphere  where  reactant  concentrations  can  be  extremely 
low,  depletion  may  be  the  limiting  factor  as  discussed  in  the  conclusions  of  this  report. 

3.  DESCRIPTION  OF  INSTRUMENTS  AND  PROCEDURES 

In  October  and  November  of  1990,  a  series  of  experiments  were  carried  out  in  Calspan’s  590  m3 
environmental  chamber  where  the  reactants  could  be  carefully  controlled.  The  purpose  of  the  experiments 
was  twofold: 

•  to  measure  the  change  in  the  aerosol  size  distribution  caused  by  the  liquid-phase  conversion  of 
S02  to  sulfate  in  cloud  droplets  during  a  cloud  cycle,  and 

•  to  measure  gas-phase  photolysis  rate  of  DMS  and  SOz  under  controlled  conditions. 

This  report  describes  only  the  cloud  cycling  experiments,  together  with  a  comparison  of  measured  and 
modeled  results. 

3.1  Calspan  Environmental  Chamber 

The  Calspan  chamber  is  a  cylinder  9.14  m  high  enclosing  a  volume  of  590  m3  (20,800  ft3).  The 
1.25-cm-thick  steel  walls  of  the  chamber  are  designed  to  withstand  pressure  differentials  of  60  kPa 
(kilopascals)  (600  mb);  however,  the  Pyrex  cover  plates  of  the  irradiation  system  currently  limit  pressure 
differentials  to  about  40  cm  of  water  (39  mb).  The  maximum  pressure  differential  was  only  about  20  mb 
in  the  set  of  experiments  discussed  here.  The  inner  chamber  surface  is  coated  with  a  fluoroepoxy  type 
urethane  developed  at  NRL,  which  has  surface  energy  and  reactivity  properties  comparable  to  those  of 
FEP  Teflon.  Illumination  for  photochemical  experiments  within  the  chamber  is  provided  by  28.6  kW 
of  fluorescent  blacklight  and  sunlamps  installed  inside  24  lighting  modules  and  arranged  in  eight  vertical 
channels  attached  to  the  wall  of  the  chamber. 

In  1990,  prior  to  this  set  of  experiments,  the  ultraviolet  irradiation  system  was  refurbished  by 
replacing  both  the  fluorescent  lights  and  black  lamps.  Subsequent  to  the  experiments,  measurements  were 
made  of  the  relative  spectral  irradiance  in  the  chamber.  The  spectrum  showed  a  peak  irradiance  at  350 
nm  and  half-power  points  at  ±20  nm. 

The  chamber  light  intensity  as  monitored  by  photodissociation  of  NOz  in  air  gave  a  k,  =  (NC^),,'1 
(d03/dt)lm0  (Wu  and  Niki  1975)  of  0.39  min1.  This  value  compares  to  kx  =  0.325  min1,  with  probable 
error  of  ±15%,  computed  from  the  photostationary  state  of  [NO],  [03],  and  [N02]. 

Absolute  aerosol  filters  are  used  for  removal  of  aerosol  particles  before  and  after  the  activated 
charcoal  filters  (300  lbs  of  charcoal).  The  air  handling  capability  made  it  possible  to  reduce  the  Aitken 
nucleus  count  to  less  than  10  cm'3  and  to  remove  those  gaseous  contaminants  which  can  be  removed  by 
charcoal  filtering.  A  tube  with  many  small  holes  around  the  top  of  the  chamber  was  used  to  direct  small 
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streams  of  distilled  water  down  the  sides  of  the  Teflon-coated  walls  to  humidify  the  chamber.  A  large 
(low  RPM)  fan  in  the  chamber  provides  circulation  (when  desired)  to  maintain  uniform  conditions  within 
the  chamber.  Figure  6  shows  the  interior  of  the  chamber.  The  chamber  has  an  airlock  for  personnel  and 
large  instrumentation  as  well  as  several  observation  and  sampling  ports.  The  temperature  of  the  chamber 
is  largely  determined  by  the  temperature  within  the  building.  The  air-handling  system  is  such  that  the 
half-time  for  reducing  aerosol  concentrations  and  trace  gas  concentrations  is  about  20  min. 


Fig.  6  —  Interior  of  Calspan’s  590  environmental  chamber 


3.2  Aerosol  Instrumentation 

The  NRL  Differential  Mobility  Analyzer  (DMA)  was  built  at  NRL  and  is  described  in  detail  in 
Hoppel  (1978)  and  Hoppel  et  al.  (1983).  The  DMA  system  has  since  been  modified  to  operate  in  a 
scanning  mode  (Wang  and  Flagan  1990)  with  an  inversion  scheme  that  provides  44  fixed  channels 
between  0.005  and  0.531  microns  radius.  A  TSI  model  3022  Condensation  Particle  Counter  (CPC)  with 
butanol  working  fluid  is  used  to  count  the  transmitted  particles.  A  HP  9845C  computer  is  used  for 
control,  data  acquisition,  and  inversion  tasks. 

The  DMA  aerosol  sample  intake  was  approximately  one  meter  from  the  chamber  wall,  and  the 
sample  immediately  passed  through  a  diffusion  dryer  to  shrink  the  particles  to  their  size  at  about  80% 
relative  humidity.  The  reduction  in  the  size  of  a  sulfuric  acid  particle  when  the  relative  humidity  is 
lowered  from  100%  to  80%  is  more  than  a  factor  of  two  in  radius  and  a  factor  of  eight  in  mass.  The 
DMA  sheath  air  was  also  dried  as  it  was  recycled,  maintaining  a  dew  point  of  approximately  4°  C. 
Sampling  from  the  chamber  at  over-pressure  with  the  DMA  required  the  CPC  makeup  air  port  and 
butanol  supply  bottle  vent  to  be  plugged  to  maintain  proper  flow  (5  cm3  s'1)  through  the  CPC  and  stop 
air  from  escaping  by  bubbling  through  the  butanol  supply  line. 
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Cloud  droplet  size  spectra  were  measured  with  a  Particle  Measuring  Systems  Inc.,  Classical 
Scattering  Aerosol  Spectrometer  Probe,  Model  100  (PMS  CSASP-100)  located  within  the  chamber 
approximately  two  meters  from  the  wall.  Droplet  spectra  were  obtained  every  20  s  using  only  the  largest 
CSASP  range  which  consisted  of  15  channels,  specified  to  be  0.6  /im  in  width  and  to  span  the  range  from 
1.0  to  10.0  /an  in  radius.  The  factory  calibration  used  two  sizes  of  glass  spheres  (3  to  5  and  10  to  15 
/an  dia.)  with  complex  index  of  refraction  m  =  1.51-0i  (at  X  =  0.5893  /an)  to  calibrate  this  range  of  the 
instrument,  which  uses  a  He-Ne  laser  operating  at  a  wavelength  of  0.6328  /an.  Glass  beads  and  water 
have  different  refractive  indices.  Based  on  the  work  of  Pinnick  and  Auvermann  (1979),  errors  of  20% 
and  larger  are  anticipated  for  certain  radii  in  the  size  range  of  interest  due  to  the  differences  in  the  index 
of  infraction. 

A  second  TSI  model  3022  CPC  monitored  the  condensation  nuclei  concentration  in  the  chamber, 
which  was  recorded  every  10  s  together  with  temperature,  pressure,  and  dew  point. 

3.3  Chamber  Equipment  and  Instrumentation 

3. 3. 1  Ozone  Generator 

The  ozone  generator  is  a  Model  T-408  ozonator  built  by  Wesbach  Ozone  System  Corp.  and  now 
handled  by  Polymetrics  of  Pleasanton,  CA.  Ozone  was  generated  in  these  experiments  by  feeding  oxygen 
rather  than  air  into  the  ozonator.  The  amount  of  ozone  generated  from  02  is  between  1.7  and  8%  by 
weight  depending  on  the  02  flow  rate,  primary  voltage,  and  temperature  of  the  cooling  medium.  Thus 
the  ozone  generator  supplies  only  03  and  residual  02  to  the  chamber. 

3. 3. 2  Sulfur  Dioxide  Monitor 

Sulfur  dioxide  was  monitored  with  a  Thermo  Environmental  Instruments  Inc.  Model  43S  High 
Sensitivity  Pulsed  Fluorescence  S02  Analyzer.  This  instrument  has  a  lower  detectable  limit  of  0.1  ppb 
S02  and  a  rejection  ratio  for  NO  of  approximately  50:1. 

3.3.3  Liquid  Water  Content  Measurement 

The  liquid  water  content  of  the  fogs  generated  in  the  chamber  was  measured  via  extinction  at  an 
infrared  wavelength  of  1 1  /im  (Chylek  1978).  The  IR  beam  enters  the  chamber  with  intensity  I,,,  traverses 
the  chamber  diameter,  is  reflected  at  the  far  wall,  and  returns  to  the  detector  (a  total  distance  of 
approximately  19  m)  where  the  intensity  I  is  measured.  The  liquid  water  content  is  then  computed  from 
the  extinction  as  determined  by  I/I0. 

3.3.4  Ion  Chromatography 

Ion  chromatography  was  performed  on  filter  samples  to  detect  sulfate,  nitrate,  ammonium  ions,  and 
methane  sulfonic  acid.  The  filters  are  90  mm  in  diameter  and  were  exposed  at  a  nominal  flow  rate  of 
1  m 3  min'1  for  10  mm,  thus  giving  a  total  volume  of  approximately  10  m 3.  The  collected  aerosol  was 
then  dissolved  from  the  filter  with  deionized  water  and  the  resulting  solution  was  injected  into  an  ion 
chromatograph. 

Calspan  chamber  instrumentation  also  included  the  following  where  the  minimum  detectable  limit 
is  included  in  parentheses: 

•  Bendix  Model  8002  Ozone  Analyzer  (1  ppb) 

•  Bendix  Model  9 101  -B  Nitrogen  Oxide  Analyzer  for  NO  and  N02  (5  ppb) 
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•  Meloy  Model  5A260  Flame  Photometric  Total  Sulfur  Gas  Analyzer  (0.5  ppb) 

•  Bendix  Model  8201  Reactive  Hydrocarbon  Analyzer  for  methane,  total  hydrocarbon  and  reactive 
(non-methane)  hydrocarbons  (5  ppb) 

3.4  Typical  Cloud-Cycle  Procedure 

The  generation  of  the  background  aerosol,  humidification  of  the  chamber,  and  subsequent  cloud 
formation  for  each  of  the  cloud  experiments  followed  a  similar  pattern.  Overnight  the  air  in  the  chamber 
was  continuously  filtered  to  remove  aerosol  and  trace  gases.  When  the  experiment  began  in  the  morning, 
the  SO2,  03  and  NOx  concentrations  were  below  the  limits  of  detection  (0.2  ppb,  2  ppb,  and  1  ppb, 
respectively).  The  charcoal  filters  did  not  remove  methane  and  were  not  very  effective  in  removing 
non-methane  hydrocarbons  (NMHC).  NMHC  concentration  was  usually,  but  not  always,  below  the 
ambient  levels.  Known  concentrations  of  03  and  SOj  were  then  introduced  into  the  chamber.  Typical 
values  were  [SOJ  =  15  ppb  and  [03]  =  100  ppb.  This  produced  a  copious  supply  of  small  particles 
which  grew  into  the  size  range  which  could  be  detected  with  the  DMA  in  about  20  min.  The  particles 
were  then  allowed  to  grow  until  the  desired  size  distribution  was  obtained.  This  usually  took  about  1  to 
2  hours.  Over  this  growth  period,  there  was  always  significant  loss  of  ozone,  presumably  as  a  result  of 
wall  losses.  To  sustain  the  aerosol  growth  rate,  additional  Oj  was  on  most  occasions  added  to  the 
chamber;  on  a  few  occasions,  additional  S02  was  also  added.  The  growth  rate  could  be  increased 
substantially  by  adding  NMHC  in  the  form  of  propene.  The  growth  occurred  at  about  the  same  rate 
whether  or  not  the  lamps  (which  simulated  the  solar  spectrum)  were  on.  Therefore  the  lamps  were  left 
off  to  prevent  heat  accumulation  in  the  chamber.  Since  the  reaction  occurred  in  the  dark  and  NMHC 
were  always  present  at  concentrations  greater  than  several  hundred  ppb,  it  was  assumed  that  the  SO2  was 
converted  to  sulfuric  acid  aerosol  by  the  Criegee  intermediate  process  (see,  for  example,  Finlayson-Pitts 
and  Pitts  1986).  The  ozone  reacts  with  the  NMHC  forming  Criegee  intermediates  which  oxidize  the  SO*. 
Figure  7  shows  the  process  of  forming  the  initial  aerosol  size  distribution.  In  this  case,  the  beginning 
concentrations  of  03  and  S02  were  120  ppb  and  15  ppb,  respectively.  The  first  size  distribution  was 
taken  20  min  after  adding  reactants  to  the  chamber,  and  the  subsequent  size  distributions  were  taken  10, 
20,  30,  50,  and  70  min  after  the  first  size  distribution. 

As  will  be  evident  later,  the  effect  of  in-cloud  conversion  is  most  readily  observed  when  the  initial 
size  distribution  is  very  steep  in  the  size  range  containing  the  smallest  particle  to  be  activated  in  the  cloud. 

After  the  desired  size  distribution  was  obtained,  the  chamber  was  filtered  to  decrease  the 
concentrations  of  the  reactants  and  aerosols.  The  chamber  air  was  typically  filtered  for  about  40  min  to 
an  hour  and  the  concentrations  were  reduced  by  about  a  factor  of  5.  This  reduction  in  the  reactant 
concentration  slowed  further  growth  of  the  particles.  The  chamber  was  then  humidified  by  passing 
distilled,  deionized  water  down  the  wall,  forming  puddles  on  the  floor.  After  about  half  an  hour,  when 
the  dewpoint  depression  was  less  than  1°  C,  the  levels  of  S02  and  03  desired  for  the  cloud  cycling 
experiment  were  added  to  the  chamber.  The  amount  of  additional  SQ2  needed  to  obtain  the  targeted 
concentration  was  always  considerably  greater  than  that  calculated.  This  was  attributed  to  loss  of  SO^ 
to  the  walls  and  to  the  water  on  the  floor  and  walls  of  the  chamber.  (Typically,  the  amount  of  SQz 
required  to  produce  5  ppb  would  have  produced  30  ppb  in  the  chamber  if  there  were  no  surface  losses.) 
The  mixing  fan  ran  continuously  in  the  chamber  except  during  the  cloud  phase.  The  fan  was  turned  off 
just  prior  to  expansion  (cloud  formation).  When  the  desired  concentrations  of  gas-phase  reactants  were 
achieved  the  chamber  was  overpressured  to  approximately  15  mb  and  held  there  for  about  5  min;  at  the 
end  of  5  min,  the  measured  dew  point  depression  was  typically  0.6°  C. 
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Fig.  7  —  Formation  of  initial  aerosol  size  distribution  prior  to  cloud 
formation.  Beginning  S02  and  03  concentrations  were  IS  ppb  and 
120  ppb,  respectively.  First  size  distribution  was  taken  20  min  after 
adding  reactants;  subsequent  size  distributions  were  taken  at  10,  20, 
30,  50,  and  70  min  after  the  first  distribution. 


Since  the  dew  point  measurement  at  these  high  humidities  has  an  accuracy  of  about  0.3°  C,  the  dew 
point  depression  was  in  the  range  0.6°  ±  0.3°  C.  The  chamber  was  then  expanded  by  exhausting  air, 
and  a  cloud  formed.  The  expansion  was  continued  until  the  pressure  was  about  5  mb  below  atmospheric 
pressure  by  using  the  chamber  exhaust  fan.  The  expansion  rate  was  varied  between  experiments  and  was 
not  constant  with  time.  Typically,  a  15  mb  expansion  was  completed  in  about  50  s.  The  pressure  was 
recorded  digitally  every  10  s,  and  was  also  recorded  on  a  strip  chart  recorder. 

After  expansion  the  cloud  was  held  for  about  4  min  to  allow  for  liquid-phase  conversion  of  SQ^ 
During  this  time,  the  cloud  droplet  size  spectrum  was  measured  with  an  optical  particle  spectrometer 
every  30  s.  After  about  3  min,  the  cloud  would  start  to  thin  at  the  top  of  the  chamber  and  near  the  walls. 
After  holding  the  cloud  for  about  4  min,  the  chamber  would  be  recompressed  with  filtered  air  and  the 
size  distribution  would  be  measured  and  compared  to  the  size  distribution  prior  to  the  cloud. 

This  cloud  cycling  process  would  then  be  repeated  on  the  same  air  using  different  values  of  the 
reactant  concentrations  and  expansion  rates. 

4.  RESULTS  OF  CLOUD-CYCLING  EXPERIMENTS 

Table  1  lists  12  cloud  cycling  experiments  in  chronological  order  and  are  designated  as  data  sets  B 
through  N.  On  November  1,  2,  and  9  (ten  data  sets),  the  in-cloud  oxidation  of  SOj  by  03  was  studied. 
In-cloud  oxidation  of  SO2  by  was  studied  on  November  10.  On  each  of  these  days,  multiple  cloud 
cycles  were  performed  on  the  same  air.  The  cycle  number  is  listed  in  the  CCY  column.  The  table  gives 
the  measured  S02  and  03  concentrations  and  indicates  when  NH3  and  H202  were  added  to  the  chamber. 
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Gas  phase  concentrations  of  NH3  and  H202  were  not  measured.  Since  both  of  these  gases  are  very 
soluble,  most  of  the  added  gas  was  undoubtedly  absorbed  by  the  water  on  the  walls  and  floor  of  the 
chamber.  The  aerosol  mass  before  and  after  the  cloud  cycle  are  those  calculated  from  the  size 
distribution  assuming  that  the  aerosol  is  sulfuric  acid  [except  for  case  J  where  the  aerosol  is  assumed  to 
be  (NH^SOJ.  Since  the  RH  in  the  DMA  is  approximately  34%,  the  sulfuric  acid  concentration  and 
density  of  the  aerosol  are  taken  to  be  50%  and  1.4  g  cm'3,  respectively.  Control  cloud-cycling 
experiments  were  carried  out  on  October  31.  These  experiments  detected  no  observable  increase  in 
aerosol  mass  when  no  reactants  were  added  to  the  chamber  or  when  SO^,  but  no  03  was  added.  There 
was  a  definite  increase  in  the  mass  of  the  size  distribution  when  both  S02  and  03  were  added  prior  to  a 
cloud  cycle  (results  not  included  in  Table  1). 

If  significant  in-cloud  conversion  occurs,  those  particles  which  were  activated  in  the  cloud  increase 
their  acidic  mass,  causing  a  separation  (local  minimum)  in  the  size  distribution.  Thus  the  critical  radius 
separates  those  particles  which  had  been  activated  in  the  cloud  from  those  which  existed  as  interstitial 
aerosol.  From  the  location  of  the  minimum,  the  activation  radius,  the  maximum  supersaturation  in  the 
cloud  and  concentration  of  cloud  droplets  can  be  inferred.  These  values  are  listed  in  the  table  for  those 
data  sets  where  detectable  separation  occurred. 

Also  listed  in  the  table  is  the  radius  range  of  the  PMS-CSASP  channel  where  the  peak  concentration 
of  droplets  occurred  at  the  time  of  maximum  liquid  water  content  (LWC).  Two  different  radius  ranges 
are  given.  One  corresponds  to  the  radii  obtained  using  the  manufacturer’s  calibration  which  was  done 
using  glass  spheres  of  known  radii.  The  other  is  an  attempt  to  correct  the  calibration  for  the  difference 
in  refractive  index  between  glass  beads  and  water  droplets  using  the  data  from  Pinnick  and  Auvermann 
(1979).  However,  our  PMS  instrument  is  not  identical  to  that  used  by  these  authors  and  the  "corrected" 
radi',  in  our  case,  may  not  represent  any  improvement.  Both  values  are  given  as  an  indication  of  the 
possible  errors  due  to  the  uncertainties  in  refractive  index  in  the  CSASP  calibration.  In  the  radius  range 
between  4  and  5  /zm,  very  little  difference  exists  between  the  factory  calibration  and  the  corrected  one. 
However  for  all  other  channels,  the  correction  was  in  a  direction  to  lower  the  indicated  radius.  Cloud 
LWCs  determined  by  IR  extinction  measurements  are  also  given.  The  absolute  concentrations  of  cloud 
droplets  given  by  the  PMS-CSASP  were  not  reliable.  The  LWC  calculated  from  the  PMS-CSASP  droplet 
spectra  was  normalized  to  the  concentration  of  droplets  inferred  from  the  DMA  and  is  given  in  Table  1 
under  the  heading  PMS/DMA.  These  values  are  always  higher  than  those  inferred  from  the  IR  extinction 
measurement  given  in  the  table.  The  IR-inferred  LWC  is  believed  to  be  the  more  accurate  because  small 
errors  in  the  PMS-measured  radius  are  magnified  in  calculating  the  volume  of  liquid  water. 

4.1  Liquid-Phase  Conversion  of  Sulfur  Dioxide  to  Sulfuric  Acid  by  Ozone 

This  section  discusses  in  detail  the  results  of  experiments  on  in-cloud  conversion  of  SOj  by  03  for 
November  1,2,  and  9  (Figs.  8  through  15).  On  each  of  these  days,  there  was  large  conversion  of  mass 
during  the  first  cloud  cycle  resulting  in  a  dramatic  change  in  the  size  distribution  as  shown  in  Figs.  8, 

1 1,  and  15.  While  the  total  mass  conversion  shown  in  Figs.  8  and  15  (data  sets  B  and  H:  6.4  and  10.3 
/xg  m  \  respectively)  was  significantly  higher  than  in  Fig.  11  (data  set  E:  3.63  /zg  m'3),  the  mass 
converted  per  cloud  droplet  was  the  largest  for  data  set  E  where  the  cloud  droplet  radii  were  the  largest. 
The  average  mass  per  droplet  converted  on  November  1 , 2,  and  9  were  4.57,  29.9,  and  6.88  femtograms, 
respectively.  This  can  be  explained  by  the  fact  that  the  cloud  droplets  formed  during  data  set  E  were 
larger.  According  to  Eq.  (12)  (and  Fig.  1),  the  growth  of  the  CCN  should  be  controlled  primarily  by 
the  cloud  droplet  radii. 

To  understand  the  change  in  the  size  distribution  upon  subsequent  cloud  cycling,  we  must  look  more 
closely  at  the  size  distribution. 
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Fig.  8  —  Changes  in  the  aerosol  size  distribution  during  the  first 
cloud  cycle  on  November  1  (data  set  B).  S02  and  0}  concentrations 
were  6.5  ppb  and  107  ppb,  respectively. 


RADIUS  {^m) 

Fig.  9  —  Changes  in  the  size  distribution  during  the  second  cloud 
cycle  of  November  1  (data  set  C).  S02  and  03  concentrations  were 
5.5  and  68  ppb,  respectively. 
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RADIUS  (jam) 

Fig.  12  —  Three  consecutive  size  distributions 
taken  between  the  first  and  second  clouds  which 
show  that  three  new  peaks  were  formed  during 
the  first  cloud  on  November  2 


0.1  1.0  10  100 
RADIUS  (urn) 


Fig.  13  —  The  cloud  droplet  size  spectrum  during  the  first  cloud  on 
November  1 .  Two  peaks  in  the  cloud  droplet  spectra  correspond  to 
the  two  peaks  in  the  aerosol  size  distribution  shown  in  Fig.  8,  taken 
after  the  first  cloud. 
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Fig.  14  —  Same  data  as  in  Fig.  11,  but  displayed 
as  cumulative  size  distributions,  together  with 
results  of  a  third  cloud  cycle.  Numbers  on  curves 
refer  to  size  distributions  taken  after  the 
numbered  cloud  cycle. 


RADIUS  (pm) 


Fig.  15  —  Size  distributions  measured  before  the 
first  cloud  cycle  and  after  the  first  and  second 
cloud  cycles  on  November  9  (data  sets  H  and  I) 
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Figure  8  shows  the  change  in  the  size  distribution  during  the  first  cloud  cycle  on  November  1  (data 
set  B)  when  both  03  (107  ppb)  and  SOj  (7  ppb)  were  present  in  the  chamber.  There  was  a  dramatic 
increase  in  the  radii  of  particles  whose  original  radii  was  greater  than  0.04  /i m.  New  peaks  in  the  size 
distribution  occurred  at  about  0.06  and  0.08  /tm  radii.  From  Fig.  1  (and  Eq.  (12))  we  would  expect 
particles  of  about  0.07  ^m  to  be  produced  by  the  evaporated  5  n m  radius  droplets.  For  the  second  cloud, 
the  expansion  was  somewhat  slower,  resulting  in  fewer  cloud  droplets  of  larger  radius.  While  the  PMS 
data  showed  that  the  droplets  were  only  marginally  larger  (with  the  peak  occurring  in  the  same  PMS 
channel  as  the  first  cycle),  the  change  in  the  size  distribution  shown  in  Fig.  9  indicates  that  particles 
greater  than  0.07  /im  were  activated  and  grew  by  a  small  amount,  adding  0.4  /ig  m'3  of  mass  to  the  size 
distribution  and  causing  a  break  in  the  size  distribution  at  0.07  /tm  to  deepen.  While  this  behavior  is 
difficult  to  detect  in  Fig.  9,  it  is  clearly  evident  when  the  difference  in  mass  before  and  after  the  cloud 
is  plotted  as  a  function  of  radius  as  shown  in  Fig.  10.  Figure  10  shows  the  change  in  mass  for  each  size 
channel  of  the  DMA  for  both  cloud  cycles.  For  the  second  cloud  cycle,  the  growth  of  particles  in  the 
0.07  to  0.09  /tm  range  removes  them  from  that  size  range  and  moves  them  to  the  0.10  to  0.13  /on  size 
range,  causing  a  modest  increase  in  mass.  The  concentration  of  particles  greater  than  about  0.2  is  less 
than  10  cm'3  so  that  errors  due  to  counting  statistics  in  channels  above  0. 15  /tm  make  the  data  above  this 
size  unreliable.  A  third  cloud  cycle  was  executed  (data  not  shown)  using  an  expansion  rate  similar  to  that 
used  for  the  first  cloud.  The  CSASP  data  indicated  that  the  cloud  droplet  radii  spectrum  was  the  same 
as  for  the  first  cloud.  No  further  conversion  of  mass  was  observed,  presumably  because  the  acidity  of 
the  cloud  droplets  due  to  the  CCN  upon  which  they  formed  was  sufficient  to  block  further  conversion. 

Figure  1 1  shows  the  differential  size  distribution  prior  to  cloud  formation  and  after  the  first  and 
second  cloud  cycles  for  November  2.  The  oscillatory  nature  of  the  size  distribution  after  the  first  cloud 
was  real  and  not  an  artifact  of  the  DMA  data  inversion  scheme.  The  maxima  and  minima  are  in  the  raw 
data.  This  behavior  occurred  during  many  of  the  data  sets,  albeit  usually  not  to  the  extent  shown  in  Fig. 
11  (see  also  Fig.  8).  Figure  12  shows  three  consecutive  size  distributions  taken  between  clouds  1  and 
2  and  demonstrates  both  the  stability  of  the  size  distribution  in  the  chamber  and  the  high  resolution  of 
the  DMA  operating  in  the  44-channel  scanning  mode. 

The  multiple  peaks  appear  to  be  associated  with  peaks  in  the  cloud  droplet  spectrum  measured  with 
the  CASAP.  Figure  13  shows  the  cloud  droplet  size  spectrum  for  the  cloud  that  produced  the  two  new 
peaks  shown  in  Fig.  8  during  the  first  cloud  cycle.  We  do  not  know  what  caused  the  multiple  peaks  in 
the  cloud  droplet  spectra;  we  know  that  they  could  not  be  identified  with  any  features  of  the  initial  aerosol 
size  distribution  on  which  the  droplets  formed.  The  model,  to  be  discussed  later,  did  not  display  any 
multiple  peaks  in  the  cloud  droplet  spectra  even  though  it  used  the  measured  expansion  rate,  included 
increase  of  solute  by  liquid  phase  conversion,  and  demonstrated  "Ostwald  ripening"  in  the  aging  cloud, 
processes  which  we  thought  might  be  responsible  for  the  multiple  peaks  in  the  size  distribution. 

During  the  second  cloud  cycle  also  shown  in  Fig.  11,  a  faster  expansion  was  employed.  This 
produced  a  higher  peak  supersaturation  which  activated  more  particles  with  smaller  radii.  Those  particles 
which  were  activated  in  the  first  cloud  evidenced  very  little  growth  during  the  second  cloud  cycle,  but 
the  radii  of  newly  activated  particles  increased  to  about  0.08  /mi.  This  behavior  can  be  explained  by 
Figure  1 .  The  first  cloud  had  a  droplet  spectrum  which  peaked  at  about  8  /t m  and  produced  evaporated 
particles  as  large  as  0. 15  /tm,  whereas  in  the  second  cloud  the  peak  in  the  droplet  spectrum  occurred  at 
only  5.5  /mi  and  produced  particles  of  about  0.08  /mi  radius.  The  larger  nuclei  produced  in  the  first 
cloud  cycle  provided  sufficient  acidity  in  the  smaller  droplets  of  the  second  cloud  so  that  no  further 
conversion  could  occur.  Figure  14  shows  the  same  data  displayed  as  cumulative  size  distributions 
together  with  the  results  of  the  third  cloud  cycle.  The  cumulative  size  distribution  enables  display  of 
multiple  size  distributions  with  less  confusion.  In  the  third  size  distribution,  the  peak  at  larger  radii  in 
Fig.  1 1  is  preserved  but  the  peak  at  0.08  /mi  is  split  by  growth  of  the  larger  particles  associated  with  that 
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peak.  This  would  indicate  that  the  maximum  supersaturation  in  the  third  cloud  was  somewhat  less  than 
that  in  the  second  cloud  and  that  the  cloud  droplets  in  the  third  cloud  were  about  20%  larger  than  in  the 
second  cloud.  This  does  not  necessarily  contradict  the  data  in  the  table  which  indicates  that  the 
PMS-CASAP  channel  containing  the  peak  concentration  of  droplets  is  the  same  for  both  data  sets.  Closer 
examination  of  the  droplet  size  distribution  shows  that  the  droplet  spectrum  was  more  sharply  peaked  at 
the  larger  sizes  during  the  third  cloud,  whereas  during  the  second  cloud  cycle,  the  spectrum  was  skewed 
to  smaller  sizes. 

The  decrease  over  time  of  the  total  aerosol  concentration  is  the  result  of 

•  diffusional  loss  of  interstitial  particles  to  the  cloud  droplets  during  the  cloud  phase  of  the  cycle, 
and 

•  "wall  losses"  during  the  period  of  several  hours  required  to  carry  out  several  cloud  cycles. 

-  •  0 

While  the  wall  losses  are  small  for  the  size  range  of  interest  here,  they  are  larger  than  expected  for  a 
totally  quiescent  chamber  because  the  mixing  fan  was  left  on  at  all  times  except  during  the  period  of 
cloud  formation.  Also,  it  should  be  pointed  out  that  throughout  the  experiment  slow  particle  growth  was 
usually  evident  in  the  hour  or  so  between  cloud  cycles.  This  was  probably  the  result  of  continued 
gas-phase  conversion  similar  to  that  used  to  produce  the  initial  size  distribution.  All  calculations  of 
in-cloud  conversion  used  the  size  distributions  measured  just  before  and  after  the  cloud  cycle.  However, 
for  purposes  of  displaying  multiple  cloud  cycles  we  used  just  one  of  a  number  of  size  distributions  taken 
between  cloud  cycles.  For  example,  in  Figs.  1 1  and  14,  the  size  distribution  used  to  illustrate  the  size 
distribution  after  the  first  cloud  is  really  the  size  distribution  measured  just  prior  to  the  second  cloud 
cycle.  The  slow  growth  during  the  hour  and  a  half  between  clouds  1  and  2  causes  the  offset  in  the  two 
size  distributions  in  the  size  range  around  0.02  to  0.03  /xm.  The  change  in  the  size  distribution  with  time 
between  clouds  can  be  judged  by  Fig.  12  which  shows  three  size  distributions  10  min  apart. 

On  November  3,  cloud  cycling  data  with  S02  and  03  in  the  chamber  were  taken  but  very  little 
conversion  was  observed.  At  the  time  this  was  puzzling,  since  we  had  obtained  dramatic  results  on  the 
two  previous  days.  However,  subsequent  analyses  of  the  data  showed  that  on  November  3  the  initial 
concentration  was  larger  (over  1000  particles  cm'3  larger  than  0.05  fim  and  the  expansion  rates  were  such 
that  cloud  droplet  concentrations  were  high  and  droplet  radii  mostly  less  than  3.5  /tm.  Figure  1  indicates 
that  under  these  conditions  conversion  rates  should  be  small.  Therefore,  the  data  from  November  3  are 
not  included  in  this  report. 

Figure  15  shows  the  results  of  the  November  9  experiment  (data  sets  H  and  I).  In  the  first  cloud 
cycle,  about  10  /xg  m3  of  material  was  formed  with  a  peak  concentration  occurring  at  about  0.1  /xm. 
This  is  at  a  somewhat  larger  size  than  we  would  have  expected  (see  Fig.  1)  for  cloud  droplets  of  the  size 
(~4.5  /xm)  indicated  by  the  PMS-CASAP.  On  the  second  cloud  cycle,  the  initial  expansion  rate  was 
slower,  with  the  total  expansion  covering  about  the  same  time  as  the  first.  This  produced  a  lower 
maximum  supersaturation  with  fewer  cloud  droplets  of  larger  radii  (~5.5  /xm).  Only  the  larger  nuclei 
in  the  previously  formed  peak  grew,  splitting  the  peak  and  forming  a  new  peak  at  about  0.14  /xm  and 
adding  about  4.5  /xg  m'3  to  the  size  distribution.  Figure  16  shows  the  mass  conversion  for  the  two  clouds 
as  a  function  of  size  (DMA  channel)  for  both  cloud  cycles.  In  the  second  cloud,  particles  larger  than 
about  0.1  /xm  were  activated  and  mass  conversion  resulted  in  growth  to  about  0.14  /xm,  causing  a  break 
in  the  size  distribution.  Again,  under  the  assumption  that  all  aerosol  material  is  sulfuric  acid,  the  size 
attained  is  larger  than  would  be  predicted  for  conversion  in  a  5.5  /xm  cloud  droplet  (see  Fig.  1). 
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Fig.  16  —  Change  in  mass  as  a  function  of  radius  during  the  two 
cloud  cycles  shown  in  Fig.  15.  6M  is  the  mass  change  in  each 
channel  of  the  DMA  plotted  at  the  midpoint  radius. 


4.2  Addition  of  NH,  to  the  Chamber 

On  several  occasions,  gaseous  NH3  was  added  to  the  chamber  before  a  third  cycle  to  determine  if 
the  absorption  of  NH3  by  the  cloud  droplet  would  increase  the  pH  of  the  droplets  and  result  in  additional 
conversion  of  S02  by  03.  After  formation  of  the  second  cloud  on  November  9,  sufficient  NH3  was  added 
to  the  chamber  to  produce  about  1  ppm  if  it  all  remained  in  the  gaseous  state.  However,  we  believe  most 
of  the  NH3  dissolved  into  the  liquid  water  on  the  walls  and  floor  of  the  chamber.  This  water  had  been 
in  contact  with  S02  and  03  during  the  previous  two  cloud  cycles  and  had  undoubtedly  become  acidic. 
On  previous  occasions  (end  of  day  on  November  2  and  3)  when  sufficient  NH3  had  been  added  to  the 
chamber  to  produce  between  100  and  500  ppb  if  all  remained  in  the  gaseous  state,  subsequent  cloud 
cycling  showed  no  significant  change  in  the  size  distribution.  This  result  remains  puzzling;  we  can  only 
assume  that  the  added  NH3  was  depleted  by  dissolution  in  the  liquid  water  before  the  cloud  cycle  was 
executed.  As  pointed  out  earlier,  when  S02  was  added  to  the  wet  chamber,  less  than  one  third  the 
calculated  concentration  appeared  in  the  chamber  as  measured  S02.  As  the  day  progressed,  it  was  easier 
to  maintain  the  SOj  concentration.  We  interpret  this  as  implying  that  the  liquid  water  progressively  came 
more  into  equilibrium  with  the  gas  phase  reactants.  However,  prior  to  every  new  cloud  cycle,  additional 
distilled  water  was  added  to  the  chamber. 

Figure  17  shows  the  size  distributions  taken  before  and  after  the  third  cloud  cycle  on  November  9 
when  NH3  was  present  with  S02  and  03.  If  we  assume  the  aerosol  (measured  in  the  DMA)  is  dry 
(NH4)2S04  with  a  density  of  1.77,  73  ng  m'3  were  formed.  We  believe  the  very  large  conversion  rate 
was  due  to  the  presence  of  the  ammonia  which  neutralized  the  aerosol.  There  was  some  question 
regarding  the  measurement  of  S02  on  this  particular  run  but  we  estimate  the  S02  concentration  to  be 
about  6  ppb.  This  presents  a  dilemma  because  if  the  S02  concentration  were  6  ppb,  and  if  it  were  all 
converted  to  (NH4)2S04,  the  maximum  yield  would  only  be  about  35  ng  m‘3.  The  most  likely  explanation 
is  that  the  (NH4)2S04  which  existed  as  a  solution  particle  at  the  high  humidity  in  the  chamber  did  not 
evaporate  to  dry  size  in  the  DMA.  Tang  (1980)  indicates  that  (NH4)2S04  aerosol  will  not  crystallize  until 
the  relative  humidity  has  decreased  to  about  35%.  For  relative  humidities  between  30  and  40%,  the 
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hydrated  aerosol  will  have  a  radius  (mass)  about  1 .3  (2.2)  times  greater  than  the  dry  particle  (Tang  1980). 
If  this  were  the  case,  the  aerosol  measured  by  the  DMA  would  be  45%  (NH4)2S04  by  weight  and  the 
solution  droplet  would  have  a  density  of  about  1.25.  The  total  aerosol  mass  added  to  the  size  distribution 
would  then  be  about  52  fig  m‘3,  of  which  23  ng  m'3  would  be  (NH4)2S04,  well  within  the  bounds  set  by 
the  available  S02.  As  explained  earlier,  the  humidity  of  the  air  sample  enters  the  inlet  tube  through  a 
diffusion  dryer  which  lowers  the  relative  humidity  from  ~  100%  to  -85%,  at  which  point  (NH^SCX, 
particles  are  undoubtedly  solution  droplets.  As  the  air  enters  the  DMA  proper,  the  humidity  of  the 
sample  is  further  decreased  because  the  humidity  of  the  sheath  air  is  about  30%.  Whether  or  not  the 
(NH4)2S04  particles  crystallize  at  some  point  within  the  DMA  is  unknown. 


Fig.  17  —  Change  in  the  size  distribution  during 
the  third  cloud  cycle  on  November  9  after  NH3 
had  been  added  to  the  chamber 


4.3  Outside  Air 

At  the  end  of  the  day  on  November  9  the  chamber  was  flushed  with  outside  air  for  an  hour  and  40 
min  before  sealing  the  chamber.  The  chamber  is  located  in  a  rural  area  of  western  New  York  State  (near 
Springville,  NY).  It  had  been  raining  all  day  and  the  air  was  believed  to  be  quite  "clean".  The  CN 
concentration  was  2000  cm 3,  and  the  measured  gas  concentrations  were  [NOx]  =  1 1  ppb,  [03]  =  3  ppb, 
and  (S02|  =  0.5+0.5  where  the  uncertainty  for  S02  is  a  result  of  the  zero  offset  uncertainty.  For  the 
first  cloud  cycle  (data  set  K)  no  reactants  were  added  to  the  chamber  and  no  change  was  observed  in  the 
size  distribution.  Before  the  second  cloud  cycle,  S02  and  03  were  added  to  the  chamber,  increasing  the 
concentrations  to  6  ppb  and  75  ppb,  respectively.  Figure  18  shows  the  size  distribution  before  and  after 
the  cloud.  The  increase  in  mass  was  about  7  fig  m'3,  more  than  doubling  the  original  aerosol  mass. 
Because  the  size  distribution  for  outside  air  was  much  less  steep  than  it  was  for  artificially  produced 
initial  size  distributions,  the  regions  of  particle  removal  are  more  graphically  portrayed.  Figure  19  shows 
the  cloud  droplet  spectrum.  The  two  peaks  observed  in  the  size  distribution  in  Fig.  18  are  ostensibly 
related  to  the  two  peaks  in  the  droplet  size  distribution. 

It  is  clear  from  Cases  B-I  and  L,  that  there  is  no  relationship  between  the  mass  conversion  in  the 
liquid  phase  and  the  liquid  water  content  for  cases  when  03  is  the  dominant  oxidant. 
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4.4  Relationship  Between  Peaks  in  the  Size  Distribution  and  Peaks  in  the  Cloud  Droplet  Spectrum 

As  shown  in  Section  2  for  liquid  phase  conversion  of  SOj  to  sulfuric  acid  by  03,  the  radial  growth 
of  the  equivalent  dry  radius  should  be  limited  primarily  by  the  size  of  the  cloud  droplet.  This  should  be 
true  for  conditions,  such  as  existed  in  our  experiment,  where  there  is  sufficient  S02  to  lower  the  pH  of 
the  cloud  droplet  (to  about  4,)  without  expending  the  S02.  Table  2  lists  the  radii  where  the  peaks 
occurred  in  the  aerosol  size  distribution  and  the  ratio  of  the  particle  radius  to  the  radius  where  the 
corresponding  peak  occurred  in  the  droplet  spectrum.  Two  sets  of  values  are  given:  one  where  the 
factory  calibration  is  used  to  determine  droplet  radii,  the  other  where  the  (uncertain)  correction  for  the 
index  of  refraction  for  water  is  used.  In  cases  where  multiple  peaks  could  be  positively  identified,  the 
values  for  the  multiple  peaks  are  also  given.  However,  it  should  be  pointed  out  that  in  every  case,  mass 
conversion  at  the  largest  radii  far  exceeded  the  mass  conversion  at  the  smaller  radii.  For  this  reason  and 
because  the  largest  peak  is  most  easily  detected  and  most  persistent  throughout  the  life  of  the  cloud,  the 
values  for  the  largest  peaks  are  believed  to  be  more  accurate.  However,  it  is  interesting  to  note  that  the 
average  ratio  for  the  second  peaks  is  very  nearly  the  same  as  the  average  ratio  for  the  largest  peaks.  The 
values  of  the  ratios  shown  in  Table  2  can  be  compared  to  the  value  (0.0137)  given  by  Eq.  (13).  While 
all  the  values  are  within  a  factor  of  two  of  each  other,  the  values  in  the  table  are  significantly  higher  than 
predicted  by  Eq.  (13).  Equation  (13)  can  be  written  as 

4  =  0.58  (k2KlK2  H0HS  P0Pst)1'9  =  0.0137,  (13') 

A  V  ' 


rp  -  Radius  where  peak  occurs  in  aerosol  size  distribution 
Rp  -  Radius  where  peak  occurs  in  cloud  droplet  spectrum 
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where  the  notation  is  the  same  as  that  given  in  Section  2.  The  last  column  of  Table  2  gives  the  factor 
by  which  the  ratio  (using  the  factory  calibration  value)  differs  from  the  predicted  ratio  of  0.0137. 
Because  of  the  ninth  root,  any  errors  in  the  product  {kJC^H^oP^^)  must  be  much  larger  than 
indicated  in  the  works  from  which  the  constants  (k2,Kl,K2,Hs,  and  H0)  were  taken  or  larger  than  expected 
errors  in  the  measurements  of  Ps,  PQ,  or  time  in  cloud.  For  example,  an  error  of  a  factor  of  20  in  the 
above  product  would  translate  to  only  a  40%  error  in  the  ratio.  Of  course,  agreement  with  the  predicted 
value  could  be  obtained  if  either  the  cloud  droplet  radii  indicated  by  the  PMS  optical  counter  were 
smaller,  or  the  aerosol  size  indicated  by  the  DMA  were  larger,  by  the  factor  given  in  the  last  column  of 
Table  2.  Possible  explanations  for  the  differences  in  the  measured  and  predicted  ratios  are  discussed 
later. 

4.5  Liquid  Phase  Conversion  of  Sulfur  Dioxide  to  Sulfuric  Acid  by  Hydrogen  Peroxide 

On  November  10,  1990,  cloud  cycling  experiments  using  H^  instead  of  03  as  the  oxidant  were 
carried  out.  Since  the  liquid-phase  oxidation  of  SO^  by  H202  is  nearly  independent  of  pH,  the  growth 
of  the  acidic  residue  (equivalent  dry  radius)  should  not  be  limited  by  the  pH  of  the  cloud  droplet  as  it  was 
in  the  case  of  the  SOj-C^  reaction.  To  inject  H202  into  the  chamber,  pure  zero  air  was  passed  over  a 
known  solution  of  H202  until  the  air  was  in  equilibrium  with  H202  and  H20  vapor  pressure  of  the 
solution  before  entering  the  chamber.  Unfortunately,  there  was  no  measurement  of  the  gas-phase 
concentration  of  H202  in  the  chamber.  Since  H202  is  very  soluble,  most  of  the  H^  is  believed  to  have 
dissolved  in  the  liquid  water  on  the  walls  and  floor  of  the  chamber.  F  the  first  cloud  (data  set  M), 
H202  was  injected  into  the  chamber  at  a  rate  which  would  produce  about  4  ppb  per  minute,  for  about  20 
min.  This  would  have  produced  80  ppb  of  H202  in  the  chamber  had  none  been  lost  to  the  walls  and  floor 
of  the  chamber.  As  stated  earlier,  we  believe  that  only  a  small  fraction  of  the  80  ppb  actually  remained 
in  the  chamber  as  gaseous  H202.  The  15  ppb  of  03  in  the  chamber  as  shown  in  Table  1 ,  was  the  remnant 
of  the  ozone  used  to  generate  the  background  aerosol  by  the  procedure  described  earlier.  Figure  20 
shows  the  size  distribution  before  and  after  the  cloud.  For  this  case  some  of  the  aerosol  particles  grew 
beyond  the  range  of  the  DMA  (greater  than  0.5  /xm)  so  that  the  mass  conversion  of  23  /xm  m 3  shown 
in  Table  1  represents  the  lower  limit  to  the  mass  conversion.  The  concentration  of  cloud  droplets  was 
considerably  lower  than  in  most  of  the  previously  described  data  sets.  The  average  conversion  per  cloud 
droplet  was  about  240  femtograms  which  is  an  order  of  magnitude  higher  than  the  highest  conversion  per 
drop  for  S02  oxidation  by  03,  and  even  about  4  times  greater  than  Case  J  where  NH3  was  used  to 
neutralize  the  acidity.  A  particle  having  a  radius  of  about  0.04  and  just  activated  in  the  cloud  increased 
in  radius  to  0.4  /xm  in  an  8  /xm  cloud  droplet.  Figure  21  shows  the  mass  conversion  as  a  function  of 
radius.  The  H202  concentration  required  for  the  observed  growth  must  be  significantly  greater  than  1 
ppb  as  shown  in  Fig.  5.  If  we  assume  that  the  H202  concentration  was  10  ppb,  then  the  ordinate  in  Fig. 
5  would  increase  by  a  factor  of  101/3  and  the  observed  increase  in  radius  would  agree  with  that  predicted 
by  Eq.  (23). 

Since  the  radii  of  the  larger  aerosol  particles  after  the  first  cloud  were  too  large  to  measure  with  the 
DMA,  an  attempt  was  made  to  remove  them  by  forming  a  cloud  while  the  mixing  fan  was  running,  thus 
impacting  the  cloud  droplets  on  the  walls  and  fan  blade.  This  effort  was  not  very  effective. 
Nevertheless,  a  second  cloud  cycle  was  attempted  and  the  results  are  given  as  data  set  N  in  Table  1 .  For 
the  second  cloud,  only  about  half  as  much  H202  was  introduced  into  the  chamber  as  in  the  prior  cloud 
cycle.  Because  of  the  large  number  of  larger  particles  existing  prior  to  the  cloud  cycle  and  the  possibility 
that  a  significant  number  of  these  particles  grew  too  large  to  be  detected  by  the  DMA,  the  interpretation 
of  data  set  N  is  uncertain. 
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Fig.  20  —  S12C  distribution  before  and  after 
cloud  cycle  when  H202  and  S02  were  reactants  in 
the  chamber  (data  set  M) 


RADIUS  (tim) 

Fig.  21  —  Mass  conversion  as  a  function  of 
radius  for  the  case  shown  in  Fig.  20.  5M  is  the 
mass  change  in  each  channel  of  the  DMA  plotted 
at  the  midpoint  radius. 
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5.  MODELING 

Since  cloud  formation  and  concentrations  of  reactants  in  the  chamber  can  be  controlled  and 
measured,  the  results  described  in  Section  4  should  be  a  good  data  set  with  which  to  compare  modeled 
results.  The  following  model  used  the  measured  expansion  rate  in  the  chamber  and  the  measured  initial 
aerosol  size  distribution  (CCN  ■spectrum)  to  predict  the  cloud  droplet  spectrum.  The  liquid-phase 
oxidation  of  S02  by  03  (or  H202)  in  the  cloud  droplets  as  giver,  by  the  equations  derived  in  Section  2 
are  used  interactively  to  predict  the  change  in  the  "dry"  aerosol  size  distribution. 

5.1  Description  of  Model 

The  change  in  the  saturation  ratio  dS  in  the  chamber  produced  by  a  given  change  in  the  pressure  dP 
is 


where  dM  is  the  mass  of  liquid  water  (per  cm'3)  formed  during  the  change  dP,  L  is  the  latent  heat  of 
water,  cp  and  R'  are  the  heat  capacity  and  gas  constant  for  air  in  cal  g  l  AT1,  p>ir  is  the  density  of  air,  pj 
is  the  water  saturation  vapor  density,  and  T  is  the  temperature.  Substituting  in  the  numerical  values  and 
taking  T  -  293  K,  the  change  in  the  saturation  ratio  is 

dS  =  -4.196  ^  -  1.841  x  105  dM.  (28) 

P 

The  first  term  gives  the  change  in  the  supersaturation  if  no  liquid  water  were  formed  (as  in  an  infinitely 
fast  expansion),  and  the  second  term  accounts  for  the  changes  in  water  vapor  concentrations  and  latent 
heat  during  condensation  or  evaporation  of  droplets.  The  expansion  is  assumed  to  be  adiabatic,  i.e.,  no 
heat  is  transferred  through  the  walls. 

The  formation  of  liquid  water  is  governed  by  the  droplet  growth  equation.  In  the  model,  the 
following  form  (Fletcher  1962)  of  the  droplet  growth  equation  is  used: 

dRi  =  DPw  s  _  A_  +  Bo  ridty  (29) 

dt  Pw  Ri  Ri  R? 

where  /?,  is  the  radius  of  a  droplet  in  the  i,h  size  class,  S  is  the  supersaturation,  D  the  diffusion  coefficient 
for  water  vapor,  ridry  is  the  dry  radius  of  the  nucleus  on  which  the  droplets  in  the  ith  class  formed,  and 
pw  is  the  density  of  water.  The  terms  containing  A  and  B0  account  for  the  effects  of  surface  tension  and 
the  vapor  pressure  lowering  due  to  the  solute  in  the  CCN.  In  our  experiments,  the  composition  of  the 
nucleus  is  well  defined  so  that  Ba  will  be  the  same  for  all  particles.  The  total  change  in  liquid  water  is 
obtained  by  summing  the  change  over  all  droplet  size  classes 

dM  =  4 t rpw  LZtRf  f^l  dt  (30> 

i  *  [dt 

where  Z,  is  the  number  of  cloud  droplets  in  the  ilh  size  class. 
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In  the  model,  the  size  bins  are  defined  by  the  midpoint  radii  of  the  44  channels  of  the  DMA;  the  bin 
boundaries  are  allowed  to  grow  (or  shrink)  as  would  a  particle  of  the  same  radius  as  the  boundary.  Thus 
a  particle  must  always  remain  in  the  same  size  bin.  The  calculation  is  started  at  saturation  (5=1)  where 
the  particle  (boundary)  radius  is  given  by 


(31) 


where  ridfy  is  the  radius  at  zero  relative  humidity  (pure  sulfuric  acid).  In  the  case  of  the  DMA  the 
measurements  were  made  at  about  34%  RH  where  the  particles  are  50%  sulfuric  acid  by  weight  and  have 
a  density  of  about  1.4.  When  the  expansion  begins,  the  radii  marking  the  channel  boundaries  grow  in 
accord  with  Eq.  (29)  and  the  change  in  liquid  water  ( dM)  is  calculated  from  Eq.  (30)  and  fed  back  into 
Eq.  (28)  at  the  next  time  step.  Values  of  the  pressure  in  the  chamber  were  recorded  every  10  s  and  these 
values  are  used  to  calculate  the  expansion  rate  required  in  Eq.(28).  Figure  22  shows  the  growth  of  the 
channel  boundaries  and  the  evolution  of  the  supersaturation  with  time  for  a  typical  expansion.  Only  the 
largest  26  of  44  channels  are  shown.  The  evolution  of  the  supersaturation  is  also  shown.  The  separation 
between  the  radii  of  unactivated  and  activated  particles  is  clearly  evident  as  the  activated  droplets  grow. 
After  the  expansion  is  completed  and  the  chamber  is  held  at  reduced  pressure,  the  total  liquid  water 
remains  constant.  The  end  of  the  expansion  is  reflected  by  a  small  discontinuity  in  the  supersaturation. 
One  channel  that  was  activated  evaporated  after  the  expansion  was  completed  about  100  s  into  the  cloud. 
As  additional  water  vapor  became  available  from  the  evaporating  droplets,  the  supersaturation  increased 
slightly.  The  liquid  water  from  the  evaporating  droplets  is  transferred  to  the  larger  droplets  (Ostwald 
ripening)  and  the  supersaturation  is  determined  by  the  supersaturation  required  to  sustain  the  smallest  size 
cloud  droplets.  Figure  23  shows  the  evolution  of  the  computed  cumulative  droplet  distribution  at  30  s, 
1,  2,  3  and  4  min  after  the  beginning  of  the  expansion  for  the  simulation  of  data  set  H.  The  measured 
size  distribution  at  34%  RH  is  shown,  as  is  the  size  distribution  adjusted  to  saturation  (dotted  line). 

The  agreement  between  the  measured  and  calculated  droplet  distributions  was  usually  quite  good  with 
respect  to  size.  However,  the  concentrations  measured  with  the  PMS  instrument  were  always  larger  than 
the  computed  value,  much  larger  than  the  concentration  required  to  obtain  agreement  with  the  total  liquid 
water  measured  by  Calspan’s  infrared  extinction  measurement,  and  larger  than  the  concentrations  inferred 
from  the  breaks  in  the  aerosol  size  distribution  as  discussed  in  Section  3.  We  therefore  accept  the 
accuracy  of  the  normalized  PMS  size  distribution,  but  not  the  absolute  value  of  concentration.  This 
would  imply  that  the  manufacturer-determined  sensing  volume  or  air  flow  of  the  PMS  was  inaccurate. 
Since  all  of  the  expansions  were  in  the  8  to  16  mb  range,  the  maximum  liquid  water  content  for  an 
adiabatic  expansion  (starting  at  saturation)  must  be  in  the  range  of  0.18  to  0.36  g  m'3  (obtained  directly 
from  Eq.  (28)  by  setting  dS  =  0).  Since  the  chamber  was  slightly  subsaturated  at  the  start  of  the 
expansion,  the  liquid  water  contents  must  be  below  these  maximum  values.  The  values  given  by  the 
Calspan  IR  liquid  water  measurement  were  always  below  the  maximum  value,  but  the  values  calculated 
from  the  PMS  droplet  spectra  normalized  to  the  concentration  implied  from  the  break  in  the  DMA  size 
distribution  exceeded  the  maximum  theoretical  value  on  several  occasions.  The  modeled  values  of  liquid 
water  were,  of  course,  always  correct  for  the  given  expansion  and  served  as  an  integral  check  on  the 
accuracy  of  the  numerical  solution  (Eq.  (28)). 
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SUPERSATURATION  (%) 

Fig.  22  —  Growth  of  channel  boundaries  during 
cloud  formation  for  a  typical  expansion.  S  is  the 
supersaturation.  This  particular  example  was  for 
data  set  B. 


RADIUS  (nm) 


Fig.  23  —  Evolution  of  cumulative  size 
distribution  as  the  larger  aerosols  are  activated 
and  grow  to  cloud  droplet  sizes.  The  initial  size 
distribution  is  the  measured  dry  size  distribution. 
Smaller  unactivated  particles  are  plotted  at  both 
their  "dry"  size  and  size  at  saturation  (S  =  0). 
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Using  the  results  of  the  previous  section  on  aerosol  growth  laws,  it  is  straightforward  to  add  the 
liquid  phase  conversion  of  S02  to  the  model.  Equation  (5),  for  oxidation  by  ozone,  or  Eq.  (20),  for 
oxidation  by  H^,  can  be  integrated  at  each  time  step  to  get  the  mass  of  sulfuric  acid  formed  within  the 
droplets.  If  vapor  depletion  is  significant,  then  Eq,  (26)  is  also  incorporated  into  the  calculation.  The 
mass  added  due  to  liquid-phase  conversion  is  used  to  recalculate  the  "dry  radius"  and  the  solute  (last) 
term  of  Eq.  (29)  at  each  time  step.  Figure  23  also  shows  the  growth  of  the  "dry"  radius  due  to 
liquid-phase  conversion  in  the  cloud  droplets  for  Case  B,  to  be  discussed  in  the  next  section. 

The  small  radii  of  the  particles  prior  to  activation  (and  the  small  radii  of  the  unactivated  particles) 
made  it  necessary  to  use  a  very  small  time  step  in  the  integration  of  the  droplet  growth  equation. 
Numerical  solutions  simulating  4  to  5  min  of  in-cloud  conversion  required  about  an  hour  to  execute  on 
a  486-type  personal  computer  running  at  33  MHz.  The  program  was  written  in  compiled  version  of 
Quick  Basic.  To  speed  up  the  computations,  aerosol  particles  which  were  several  channels  smaller  than 
the  smallest  size  for  activation  were  ignored  and  they  remained  at  their  size  at  saturation.  It  was  verified 
that  this  assumption  had  no  effect  on  the  solution;  i.e.,  the  amount  of  additional  water  taken  up  by  the 
unactivated  particles  during  the  cloud  cycle  is  negligible,  as  is  the  liquid-phase  conversion  in  unactivated 
particles. 

5.2  Model  Results  for  Specific  Runs 

The  results  of  the  modeling  are  summarized  in  the  second  part  of  Table  1.  For  comparing  how  well 
the  model  simulated  the  measured  parameters,  the  important  parameters  are  the  maximum  supersaturation, 
and  size  and  concentration  of  cloud  droplets.  As  indicated  earlier,  the  change  in  the  size  distribution  can 
be  used  to  indicate  the  maximum  supersaturation  in  the  cloud  for  those  cases  where  measurable  in-cloud 
conversion  occurred.  These  experimentally  determined  supersaturations  (Table  1(a))  can  be  compared 
to  the  model  predicted  supersaturations  (Table  1(b)).  The  modeled  and  experimentally  determined  values 
of  the  maximum  supersaturation  track  quite  well  but  the  modeled  values  are  always  somewhat  higher. 
There  is  relatively  good  agreement  in  the  measured  and  modeled  cloud  droplet  size  and  this  is  reflected 
in  Table  1  where  the  peaks  in  the  modeled  and  measured  spectrum  correlate  quite  well.  Figure  24  shows 
a  typical  comparison  of  measured  and  modeled  cloud  droplet  size  spectra  at  2,  3,  4,  and  5  min  into  the 
cloud  (for  data  set  B).  Six  measured  droplet  spectra  were  taken  during  the  same  time  period  but  some 
are  indistinguishable  from  each  other  (there  appear  to  be  only  two)  in  Fig.  24.  Ostwald  ripening  is 
clearly  evident  in  the  modeled  result  where  growth  of  the  larger  particles  continued  even  though  the  liquid 
water  content  is  constant.  The  agreement  is  better  using  the  radii  from  the  factory  calibration  of  the 
PMS-OPC  rather  than  the  radii  from  the  attempt  to  correct  for  the  refractive  index  of  water.  The  size 
and  number  of  droplets  determines  the  total  cloud  liquid  water  content.  As  stated  earlier,  the 
concentration  of  droplets  measured  with  the  PMS-OPC  is  too  large;  only  the  radial  distribution  is  assumed 
to  be  valid.  The  droplet  radii  intervals  listed  in  Table  1(b)  are  the  boundary  channels  containing  the 
maximum  number  of  cloud  droplets  at  the  time  when  the  cloud  has  its  maximum  liquid  water  content, 
which  is  also  given  in  Table  1(b). 

We  believe  the  agreement  between  the  computed  and  measured  cloud  parameters  is  as  good  as  can 
be  expected  considering  the  measurement  errors  involved  in  the  PMS-OPC  and  considering  the  difficulty 
in  selecting  the  critical  radius  by  examining  the  change  in  the  size  distribution  before  and  after  cloud 
conversion.  Any  errors  in  selecting  the  size  of  the  smallest  particle  activated  affects  the  experimentally 
determined  cloud  droplet  concentration  and  supersaturation  given  in  Table  1(a).  The  total  number  of 
cloud  droplets  as  determined  from  the  break  in  the  measured  size  distribution  and  the  number  predicted 
by  the  model  show  significant  differences.  This  can  be  attributed  to  the  steepness  of  the  initial  size 
distributions.  As  a  result,  the  concentration  of  particles  activated  in  the  model  is  very  sensitive  to  the 
expansion  rate  and  to  the  amount  of  acid  in  the  particles  at  the  reference  humidity.  Even  though  our 
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Fig.  24  —  Comparison  of  measured  and  modeled  cloud  droplet 
spectra.  Modeled  results  are  plotted  for  2,  3, 4,  and  5  min.  Measured 
spectra  are  six  spectra  taken  during  the  same  time  period  and  several 
lie  on  each  other  making  them  impossible  to  resolve  individually. 
Example  is  from  simulation  of  data  set  B. 


DMA  was  operating  at  very  high  resolution  (44  channels),  the  steepness  of  the  size  distribution  would 
often  result  in  a  factor  of  two  change  in  the  total  cloud  droplet  population  if  one  additional  channel  were 
activated. 

The  largest  errors  in  the  modeled  cloud  parameters  stem  from  our  inadequate  knowledge  of  the  exact 
relative  humidity  in  the  chamber  prior  to  the  expansion  (due  to  finite  accuracy  of  the  dew  point 
hygrometer)  and  to  the  rather  coarse  data  on  the  expansion  rate  (pressure  measurement  every  10  s). 
Since  the  chamber  is  slightly  subsaturated  (by  an  unknown  amount),  a  curve-fit  to  the  pressure  points  was 
used  to  extrapolate  back  to  the  starting  point  of  the  expansion.  In  our  calculations,  it  is  assumed  that 
saturation  was  achieved  after  5  s  of  expansion,  at  which  time  model  calculations  began  assuming  a 
saturated  chamber.  As  stated  earlier,  the  duration  of  the  expansion  varied,  but  typically  lasted  about  1 
to  1.5  min. 

5.2. 1  Modeled  Results  for  Oxidation  of  Sulfur  Dioxide  by  Ozone 

The  concentrations  of  reactants  were  measured  only  for  the  cases  in  which  liquid-phase  oxidation 
of  SOs  was  by  03.  Therefore,  only  for  data  sets  B,  C,  E,  F,  G,  H,  I,  and  L,  can  a  comparison  of 
modeled  and  experimental  results  be  meaningful.  As  was  discussed  in  the  last  section,  the  most  dramatic 
changes  occurred  on  the  first  cloud  cycle  on  each  of  the  three  days  when  SO*  oxidation  by  03  was 
studied.  The  results  of  all  these  runs  are  given  in  Table  1(b),  but  detailed  results  will  be  shown  only  for 
data  sets  B,  E,  and  H,  which  were  the  first  cloud  cycles  of  each  day,  and  data  set  L,  which  is  data  in  air 
which  contained  natural  aerosols  as  the  starting  size  distribution. 

The  modeled  change  in  the  dry  size  distribution  after  30  s,  1,  2,  3,  4,  and  5  min  are  shown  in  Figs. 
25  through  28  for  data  sets  B,  E,  H,  and  L.  The  double-peaked  characteristic  caused  by  the  mass 
conversion  in  cloud  droplets  is  clearly  present,  but  the  total  mass  conversion  is  much  smaller  than  the 
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observed  mass  conversion  shown  in  Figs.  8,  11,  15,  and  18.  For  case  E  (Fig.  26),  particle  growth  was 
largest  because  the  cloud  droplets  were  considerably  larger  than  in  cases  B  and  H.  For  data  set  L  (Fig. 
28),  the  case  when  the  aerosol  was  ambient  aerosol  from  flushing  the  chamber  with  outside  air,  the  initial 
aerosol  was  assumed  to  be  nonacidic,  but  to  have  water  absorption  characteristics  similar  to  sulfuric  acid 
aerosol.  In  all  other  cases,  the  initial  aerosol  was  assumed  to  be  acidic. 

Figure  29  shows  the  cumulative  (dry)  size  distributions  at  0,  30  s,  1,  2,  3,  and  4  min  for  data  set 
H.  This  is  the  same  dry  data  as  shown  in  Fig.  27  as  a  differential  size  distribution.  Also  shown  in  Fig. 
29  is  the  cumulative  measured  size  distribution  after  the  cloud  showing  the  large  discrepancy  between  the 
measured  and  modeled  conversion. 


Fig.  29  —  Cumulative  distribution  for  simulation 
of  data  set  H  including  the  measured  size 
distribution  after  the  cloud 


It  is  clear  from  Table  1  that  the  modeled  change  in  the  mass  during  the  cloud  cycles  for  these  cases 
(S02  oxidation  by  03)  was  much  lower  than  that  actually  measured;  sometimes  the  modeled  mass 
conversion  was  more  than  an  order  of  magnitude  smaller  than  observed.  This  is  consistent  with  the  results 
shown  in  Table  2  and  the  discussion  in  Section  4.4.  When  the  constant  in  the  second  term  of  Eq.  (5)  was 
increased  sufficiently  to  make  the  ratio  of  r/R  given  by  Eq.  (13')  to  agree  with  the  value  of  the  ratio 
given  in  Table  2  (instead  of  the  value  0.0137  given  in  Eq.  (13)),  then  the  agreement  between  the  modeled 
mass  conversion  and  observed  mass  conversion  were  quite  good,  as  would  be  expected.  The  added  mass 
using  the  adjusted  constant  is  shown  in  the  thirteenth  column  of  Table  1(b)  for  data  sets  B,  E,  and  H. 
However,  this  "forcing"  of  agreement  between  observation  and  modeling  is  unsatisfactory,  since,  as 
discussed  in  Section  4.4,  the  magnitude  of  the  required  change  in  the  constant  is  larger  than  any 
anticipated  errors  in  the  experimentally  determined  parameters  which  make  up  this  constant  (in  Eq.  (5)). 
It  does,  however,  show  that  the  inferences  of  our  simpler  arguments  given  in  Section  2.1,  based  only  on 
the  growth  law  and  cloud  droplet  size,  agree  with  the  more  rigorous  model  calculations. 

The  reason  for  the  large  discrepancy  in  mass  conversion  between  the  measured  and  modeled  results 
is  not  known.  It  may  be  that  the  ozone  generator  (described  in  Section  3.3)  which  uses  a  high-voltage 
discharge,  may  have  produced  small  concentrations  of  some  unknown  product  such  as  H202  or  some 
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excited  molecule  which,  when  dissolved  in  the  cloud  droplets,  produced  spurious  oxidation  of  S02.  For 
example,  if  in  addition  to  producing  82  ppb  of  03  the  generator  also  produced  2  ppb  of  H202  in  the 
chamber,  then  the  resulting  conversion  of  S02  would  have  produced  10.4  ^g  m'3  instead  oi  0  976  /xg  m'3 
of  sulfuric  acid  aerosol  for  the  case  of  data  set  H.  The  change  in  the  modeled  size  distribution  for  data 
set  H  assuming  2  ppb  of  H202  is  shown  in  Fig.  30  where  it  is  compared  to  the  measured  change  in  the 
size  distribution.  This  result  is  obtained  by  adding  the  conversion  rate  equations  for  03  and  H202 
oxidation  of  S02  (Eqs.  (5)  and  (20))  in  the  model.  Even  though  the  concentration  cf  Qj  is  82  ppb  and 
H202  is  only  2  ppb,  the  conversion  is  due  primarily  to  the  H202  and  is  limited  by  the  availability  of 
H202.  This  calculation  (Fig.  30)  is  given  to  show  that  small  amounts  of  unknown  oxidants  could  produce 
the  observed  conversion.  There  is  no  way  to  know  for  certain  what,  if  any,  trace  oxidants  are  formed 
by  the  O,  generator  without  extensive  tests  on  the  output  of  the  generator  und^r  the  conditions  of  our 
experiments,  where  the  output  was  introduced  directly  into  a  saturated  chamber.  The:  *'  ‘ests  are  beyond 
the  scope  of  this  investigation. 


Fig.  30  —  Modeled  growth  in  the  dry  size 
distribution  for  iLta  set  H  when  2  ppb  of  H202  is 
ass  :  ned  to  be  present  in  addition  to  82  ppb  of 
03.  Measured  size  distribution  after  the  cloud  is 
shown  for  comparison. 


The  above  explanation  loses  some  of  its  appeal  when  it  is  remembered  that  less  conversion  occurred 
on  subsequent  cloud  cycles.  This  was  explained  in  terms  of  the  increased  droplet  acidity  due  to  the 
increased  size  of  the  CCN  on  which  the  cloud  droplet  formed  on  subsequent  cloud  cycles.  Before 
repeating  a  cloud  cycle,  the  ozone  generator  was  used  to  increase  the  03  concentration  to  the  70  to  100 
ppb  range.  This  "topping-off  of  the  ozone  concentration  required  less  ozone  to  be  injected  into  the 
chamber  than  the  initial  injection,  but  any  extraneous  product  produced  by  the  generator  would  have  also 
been  injected  before  subsequent  cloud  cycles,  and  we  would  not  have  expected  to  see  the  pH  limitation 
described  in  Section  3. 

In  a  recent  paper,  Hansen  et  al.  (1991)  suggests  that  their  data  indicate  that  oxidation  of  S02  to 
sulfate  in  a  dispersed  system  is  different  than  in  the  bulk  system,  and  that  reaction  rates  in  droplets  may 
be  quite  different  from  those  predicted  by  extrapolation  of  bulk  solution  results.  Our  results  could  be 
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interpreted  as  further  experimental  evidence  to  support  this  hypothesis.  Until  this  hypothesis  finds  better 
theoretical  foundation,  we  are  reluctant  to  invoke  it  to  explain  our  data. 

A  possible,  but  less  probable,  explanation  of  the  discrepancy  between  modeled  and  measured 
conversion  rates  is  that  some  ammonia  may  have  been  in  the  chamber.  Ammonia  used  in  the  chamber 
on  previous  runs  and  removed  by  the  charcoal  filters  could  conceivably  have  been  released  during 
filtering  to  lower  the  concentrations  of  aerosols  after  the  formation  of  the  initial  size  distribution.  No 
additional  filtering  was  done  over  the  number  of  hours  during  which  the  chamber  was  humidified  and 
several  cloud  cycles  were  executed.  Also,  for  the  outside  air  sample  (data  set  L)  the  air  was  not  passed 
through  the  charcoal  filters,  and  the  observed  conversion  was  still  much  higher  than  the  modeled 
conversion. 

A  fourth  possible  explanation  would  be  that  the  cloud  droplet  radii  were  larger  than  those  indicated 
by  the  measurement  with  the  PMS-OPC  by  the  difference  factor  in  Table  2  (factor  of  10%  to  100%). 
However,  the  rather  good  agreement  between  the  measured  and  modeled  values  of  cloud  droplet  radii  and 
the  constraint  on  the  total  liquid  water  content  make  this  latter  explanation  improbable. 

5.2.2  Modeled  Results  for  the  Case  when  Ammonia  was  Present 

Data  set  J  showed  large  conversion  of  S02  by  03  when  gaseous  NH3  was  present.  The  NH3 
presumably  dissolves  in  the  droplets  where  it  neutralizes  the  H2S04  forming  (NH4)2S04.  The  ammonium 
sulfate  keeps  the  pH  from  dropping  and  facilitates  oxidation  of  S02  by  03  according  to  Eq.  (2).  If  it  is 
assumed  that  the  presence  of  NH3  only  effects  the  pH  ([//+])  of  Eq.  (1),  then  the  conversion  rate  of  S02 
to  H2S04  to  (NH4)2S04  can  be  readily  calculated  by  the  model.  The  conversion  rate  is  then  given  by 

dM„  . 

- 2  =  K  R3  PnP v  (32) 

dt  °  * 

where  A"  is  1.1  x  1012,  1.1  x  1010,  or  1.1  x  108  for  pH  of  7,  6,  or  5,  respectively.  The  model 
calculations  were  carried  out  for  pH  of  5,  6,  and  7,  and  gave  mass  conversions  of  0.86,  32,  and  35  ng 
m  3,  respectively.  If  the  water  associated  with  the  ammonium  sulfate  aerosol  did  not  completely 
evaporate  (as  discussed  in  Section  4.2),  the  mass  conversion  as  inferred  from  the  DMA  size  distribution 
would  be  1.92,  71  and  79  /ig  m'3,  respectively.  For  the  case  when  pH  =  6  and  higher,  the  depletion  of 
S02  was  the  dominant  factor  limiting  conversion.  For  the  case  of  pH  =  7,  the  SC^  was  used  up  in  about 
10  s  (invalidating  the  assumption  of  gas-liquid  equilibrium).  For  the  case  of  pH  =  6,  about  90%  of  the 
S02  was  used  up  4  min  into  the  cloud,  and  for  pH  =  5,  2%  of  the  SO2  was  used  up  4  min  into  the  cloud. 
There  is  a  sharp  cutoff  in  the  conversion  rate  as  the  pH  drops  from  6  to  5.  The  modeled  results  for  the 
case  when  the  pH  was  assumed  to  be  6  is  shown  in  Fig.  31,  where  the  aerosol  is  assumed  to  be  dry 
(NH4)2S04  with  a  density  of  1.77.  Also  shown  is  the  measured  size  distribution  after  the  cloud. 

5.2.3  Model  Results  for  the  Oxidation  of  Sulfur  Dioxide  by  Hydrogen  Peroxide 

Liquid  phase  oxidation  of  S02  by  H202  in  a  cloud  was  attempted  on  November  11  and  the 
experimental  results  (data  set  M)  are  shown  in  Fig.  32.  As  discussed  in  Section  4.5,  no  measurements 
of  H202  concentrations  in  the  chamber  were  available,  but  H202  vapor  could  be  injected  into  the  chamber 
by  flowing  pure  air  over  a  solution  of  H202  at  a  rate  which  would  result  in  a  4  ppb-per-minute  increase. 
This  was  done  for  about  20  min  prior  to  cloud  formation  for  data  set  M.  As  stated  earlier,  we  believe 
that  only  a  small  fraction  of  the  H202  remains  in  the  vapor  state  and  most  of  the  H202  was  lost  to  the 
water  on  the  walls  and  floor  of  the  chamber  prior  to  cloud  formation.  The  model  was  run  using  Eq.  (20) 
for  the  liquid  phase  conversion  of  S02  by  H202,  3.5  ppb  of  S02,  as  measured  just  prior  to  cloud 
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Fig.  31  —  Modeled  growth  of  the  dry  size 
distribution  when  the  presence  of  NHj  was 
assumed  to  hold  the  pH  of  the  droplets  at  6.  The 
resulting  aerosol  is  assumed  to  be  dry 
(NH4)2S04.  Measured  size  distribution  after  the 
cloud  is  shown  for  comparison. 


RADIUS  (um) 

Fig.  32  —  Modeled  growth  of  dry  size  distribution  using  initial 
distribution  from  data  set  M.  The  S02  concentration  is  3.5  and  the 
concentration  of  H202  was  assumed  to  be  10  ppb.  Also  shown  is  the 
measured  size  distribution  after  the  cloud  cycle. 
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formation,  and  for  several  assumed  concentrations  of  H202.  For  10  ppb  and  2  ppb  of  H202,  the  mass 
conversions  were  26.9  and  7.42  pg  m'\  respectively.  The  case  for  10  ppb  of  H202  is  shown  in  Fig.  31 
and  the  result  given  in  Table  1(b),  case  M.  While  the  exact  concentration  of  H202  is  not  known,  it  is 
clear  that  reasonable  agreement  between  the  measured  and  modeled  size  distributions  can  be  achieved  with 
reasonable  values  of  H202.  For  the  case  of  10  ppb  of  H202,  about  2.5  ppb  of  the  H202  and  S02  were 
consumed,  whereas  in  the  case  of  2  ppb  of  H202,  about  0.73  ppb  of  the  H202  and  S02  were  consumed. 

There  appeared  to  be  increased  growth  of  the  equivalent  dry  radius  of  the  aerosol  after  S02  and  H202 
were  added  to  the  wet  (but  subsaturated)  chamber.  This  phenomenon  was  not  investigated  in  detail  and 
usually  the  cloud  was  formed  as  soon  after  H202  injection  as  possible,  so  that  the  H202  would  not  be  lost 
to  the  wet  walls  of  the  chamber.  The  growth  of  H2S04  aerosol  by  oxidation  of  S02  by  H202  in  the  high 
humidity  (subsaturated)  environment  can  be  calculated  from  Eq.  (20)  where  R  is  now  interpreted  as  the 
equilibrium  radius  of  the  solution  droplet  at  high  humidity.  The  relationship  between  the  radius  of  the 
solution  droplet  and  the  radius  ra  of  the  pure  acid  particle  after  all  water  is  removed  can  be  written  as 


S  =  Exp 


A 

R 


(33) 


provided  the  saturation  ratio  is  high  enough  (0.85  <  S  <  1)  that  the  mass  of  the  water  is  much  greater 
than  the  mass  of  the  pure  acid.  Expressing  Ma  in  Eq.  (20)  in  terms  of  the  radius  ra,  and  using  Eq.  (33) 
to  eliminate  ra  gives  a  differential  equation  for  R  which  has  the  solution 


*  =  Ci 


~  -  -  Ct  >"  (R/Ro) 


(34) 


where  Ra  is  the  radius  of  the  solution  droplet  at  t  =  0  and  at  a  supersaturation  S,  and 

-l 


C,  =  — 

B„ 


3.9  x  1010  PSPH 
4/3  7 r  pa 


(35) 


C,  =  3(ln  S) 


3.9  x  1010  PSPH 


4/3  7T  pa 

At  any  time  t,  the  wet  radius  can  be  converted  to  the  radius  ra  of  a  pure  acid  particle  by  using  Eq.  (33). 


Solutions  to  Eq.  (34)  for  saturation  ratios  of  1,  0.995,  0.99,  0.98,  and  0.95  are  shown  in  Fig.  33 
in  terms  of  the  pure  acid  radius,  for  initial  radii  of  0. 1  and  0.2  pm  and  when  Ps  =  10  ppb  and  PH  =  10 
ppb.  Significant  growth  occurs  over  a  12  hour  period  only  if  the  relative  humidity  is  above  about  98%, 
and  for  S  =  1  the  growth  is  very  rapid. 

The  pH  can  be  calculated  by  combining  Eqs.  (3)  and  (33) 

[^♦1  _  (2-436)  4/3  *  Pa 

1  J  =  ra 


R 


-  In  S 


(36) 
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TIME  (hr) 


Fig.  33  —  Growth  of  pure  acid  particle  (water  removed)  at  saturation 
ratios  of  1.0,  0.995,  0.99,  0.98,  and  0.95  in  an  environment  where 
there  is  10  ppb  of  both  S02  and  H202.  Initial  dry  sizes  are  0.1  and 

0.2  pm. 


The  pHs  calculated  from  Eq.  (36)  are  below  2,  which  was  stated  earlier  to  be  below  the  range  of  validity 
for  Eq.  (20).  However,  the  conversion  rates  fall  off  rather  slowly  with  decreasing  pH  (Martin  1984)  and 
Eq.  (20)  can  be  considered  correct  within  a  factor  of  2  for  1  <  pH  <2  and  within  a  factor  of  5  for  0< 
pH  <1. 

This  latter  discussion  has  been  included  because  questions  regarding  the  possible  importance  of 
gas-to-particle  conversion  in  wet  aerosols  (hazes)  when  the  environmental  humidity  is  high  are  often 
raised.  Oxidation  of  SO2  by  H^  has  sometimes  been  suggested  as  an  important  mechanism  in  the 
atmosphere  because  it  can  occur  at  low  pH.  Although  high  humidities  can  occur  at  the  top  of  the 
planetary  boundary  layer,  even  when  there  are  no  visible  clouds,  and  at  ground  level  during  the  night  due 
to  radiative  cooling,  it  seems  unlikely  that  the  very  high  humidities  indicated  here  would  prevail  for  long 
enough  periods  of  time  to  make  conversion  by  H202  important  in  the  cloud-free  atmosphere. 

6.  CONCLUSIONS 

The  cloud  cycling  experiments  and  modeling  described  in  this  report  demonstrate  that  liquid-phase 
conversion  of  S02  in  cloud  droplets  can  have  a  dramatic  effect  on  the  aerosol  size  distribution. 
Observations  of  the  aerosol  size  distribution  before  and  after  the  aerosol  has  been  processed  through  a 
(nonprecipitating)  cloud  is  a  very  effective  way  of  studying  the  liquid-phase  conversion  of  trace  gases  in 
clouds.  For  the  case  when  S02  and  O3  were  introduced  into  the  chamber,  the  observed  conversion  was 
clearly  iimited  by  the  cloud  droplet  acidity  produced  by  liquid-phase  conversion  of  SOj  to  H2S04  and  by 
the  acidity  of  the  CCN  on  which  the  cloud  droplet  formed.  The  amount  of  conversion  is  a  strong 
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function  of  the  droplet  radius.  The  dilution  of  acid  in  a  large  droplet  raises  the  pH,  resulting  in  rapid 
conversion  of  S02  and  causing  increased  mass  of  the  CCN  after  cloud  evaporation. 

Starting  with  well-established  formulae  for  the  oxidation  of  S02  by  either  03  or  H202  in  Henry’s 
Law  equilibrium  with  bulk  water,  growth  laws  for  the  increase  in  "dry"  radius  were  derived.  These 
chemical  conversion  formulations  were  used  with  a  microphysical  model  of  cloud  formation  to  predict 
the  change  in  the  aerosol  size  distribution.  The  observed  conversion  for  the  case  of  SOj  oxidation  by  03 
was  considerably  larger  than  that  calculated.  The  reason  for  the  discrepancy  between  the  observed  and 
modelled  conversion  is  not  known  but  several  possibilities  were  discussed  in  Section  5.2.1.  The  cause 
of  the  difference  between  the  observed  and  predicted  conversion  rates  should  be  explored  further  in  future 
cloud  chamber  experiments. 

The  presence  of  (gaseous)  NH3,  together  with  S02  and  03,  increased  the  liquid-phase  conversion 
considerably.  This  undoubtedly  is  due  to  the  fact  that  dissolved  NH3  reacted  with  H2S04,  increasing  the 
pH  of  the  droplet.  The  increase  in  pH  promotes  SC^  absorption  and  oxidation  of  S02  within  the  droplet. 

For  the  case  of  SO2  oxidation  by  H202,  no  pH  limitation  to  the  conversion  of  S02  to  H2S04  was 
observed.  The  equivalent  "dry"  size  of  the  aerosol  increased  dramatically  and  appeared  to  be  limited  only 
by  availability  of  the  gas-phase  reactants. 

An  annoying  aspect  of  the  data  was  the  multiple  peaks  in  the  size  distribution  often  generated  in  a 
single  cloud  cycle.  While  these  seemed  to  correlate  with  corresponding  peaks  in  the  cloud  droplet 
spectrum,  the  accuracy  and  resolution  of  the  cloud  droplet  spectrum  was  not  adequate  to  verify  that 
hypothesis  conclusively.  Furthermore,  even  if  the  multiple  peaks  in  the  cloud  spectra  are  the  cause,  the 
source  of  the  multiple  peaks  in  the  cloud  spectra  is  not  understood.  Possible  sources  of  the  multiple 
peaks  could  be  related  to  the  dynamics  of  the  expansion  or  preexisting  "ripples"  in  the  size  distribution 
prior  to  cloud  formation  which  were  too  small  to  be  resolved  by  the  DMA. 

The  equations  given  in  this  report  can  be  used  to  make  some  general  comments  on  the  expected 
behavior  of  S02  concentrations  in  the  remote  marine  boundary  layer  (MBL).  In  the  following  scenarios, 
we  will  assume  that  the  flux  of  DMS  (dimethyl  sulfide)  from  the  ocean  is  140  /xg  S  m'2  day'1  (Bates  et 
al.  1987),  that  the  DMS  is  converted  uniformly  in  a  1  km  MBL  to  give  a  S02  volume  source  rate  of  3.2 
x  10-18  g  cm'3  s'1  and  that  DMS  is  the  only  source  of  S02.  It  is  also  assumed  that  the  loss  of  S02  by 
in-cloud  conversion  is  the  primary  loss  mechanism.  Ignoring  loss  by  gas-phase  conversion  and  surface 
deposition  is  obviously  an  oversimplification,  but  will  give  order-of-magnitude  estimates  for  the  effect 
of  in-cloud  conversion  on  S02  concentrations. 

Scenario  1 .  Liquid-phase  oxidation  of  S02  by  03.  Since  the  concentration  of  03  ( -  20  ppb)  is  much 
larger  than  that  of  S02,  the  03  concentration  can  be  considered  constant.  Equations  (5)  and  (26)  can  be 
written  as 


and 


dMa 

dt 


1.31  x  10'3  R9  P0  Ps 


Mi 


(37) 


dPt  dM. 

_ i  =  -228Z  _ ? 

dt  dt 


(38) 
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where  Ss  is  the  source  of  S02  from  DMS  (1 . 1  x  10'15  atm  s'1).  If  we  assume  that  the  clouds  in  the  MBL 
consist  of  100  droplets  cm'3  of  10  jim  radius  (liquid  water  content  of  .42  g  m'3),  and  that  a  parcel  of  air 
spends  one-twentieth  of  its  time  in  a  cloud,  then  the  simultaneous  solution  to  Eqs.  (37)  and  (38)  will  give 
the  S02  concentrations  and  the  aerosol  H2S04  loading.  Starting  from  zero  concentrations,  Fig.  34  shows 
the  buildup  of  S02  and  H2S04  over  the  four  days.  During  the  first  half-day,  more  of  the  SOj  is 
converted  to  H2S04  than  remains  as  S02;  however,  at  longer  times,  the  conversion  rate  slows  down  and 
more  sulfur  remains  as  S02.  This  scenario  produces  only  sulfuric  acid  aerosol  and  the  concentration  of 
S02  is  greater  after  four  days  than  is  found  in  the  remote  MBL. 


Fig.  34  —  Concentrations  of  S02  and  acidic  mass  assuming  DMS  is 
oxidized  to  S02  and  S02  is  converted  to  sulfuric  acid  by  liquid-phase 
oxidation  by  03  within  clouds.  It  is  assumed  that  the  air  spends 
one-twentieth  of  its  time  in  MBL  clouds. 


Scenario  2.  S02  is  oxidized  by  03  but  NH3  is  present  to  partially  neutralize  the  sulftiric  acid.  This 
is  a  distinct  possibility  since  typical  concentrations  of  NH3  (0.2  ppb,  Warneck  1986)  found  in  the  MBL 
are  usually  greater  than  typical  S02  concentrations  and  NH3  is  more  soluble  than  S02.  As  seen  earlier 
(Section  5.2.2),  the  conversion  rate  depends  strongly  on  the  pH  of  the  droplet.  For  fixed  pH,  the 
conversion  rate  is  given  by  Eq.  (32)  where  K  depends  strongly  on  the  pH  of  the  droplet.  Equations  (32) 
and  (38)  have  a  solution  of  the  form 


Ps  = 


~C 


(39) 


where  C  =  195  ZXF?P0  is  the  reciprocal  of  the  time  for  achieving  equilibrium.  For  pH  of  6,  5,  and  4, 
the  corresponding  time  constants  are  1 .2  hours,  5.4  days  and  540  days;  and  the  corresponding  equilibrium 
concentrations  of  SC^  are  0.005  ppb,  0.5  ppb,  and  50  ppb.  These  calculations  are  based  on  the  same 
cloud  properties  as  Scenario  1  and  a  parcel  of  air  spending  one-twentieth  of  its  time  in  a  nonprecipitating 
cloud.  It  is  clear  that  if  NH3  concentrations  are  large  enough  to  keep  the  pH  of  the  droplet  as  high  as 
6,  then  concentrations  of  S02  will  be  very  low  and  most  S02  will  be  converted  to  sulfate.  In  the  other 
extreme,  if  the  pH  is  as  low  as  4,  the  S02  concentrations  will  slowly  build  up  to  unrealistically  high 
values. 


Scenario  3.  If  oxidation  of  S02  in  the  remote  MBL  is  by  H202,  then  the  depletion  of  H202  as  well 
as  S02  must  be  included.  We  were  unable  to  identify  any  published  measurements  of  H202  in  the  remote 
MBL;  however,  measurements  at  a  number  of  locations,  mostly  over  North  America  and  Europe,  are 
tabulated  in  review  articles  by  Sakugawa  et  al.  (1990)  and  by  Gunz  and  Hoffmann  (1990).  Typical  H202 
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concentrations  in  the  MBL  are  probably  of  the  same  order  of  magnitude  as  are  SO2  concentrations,  with 
both  being  highly  variable.  In  the  low  NOx  environment  of  the  remote  MBL,  H202  may  well  be  greater 
than  S02  concentrations. 

Equations  (20)  and  (38),  which  account  for  aerosol  mass  and  S02  concentrations,  can  be  augmented 
by  the  following  equation  to  account  for  hydrogen  peroxide 


dPh 

dt 


dMa 

-228  Z - f  +  S 


dt 


H 


(40) 


where  SH  is  the  source  strength  of  H202.  Solutions  to  these  three  equations  for  the  same  cloud  properties 
as  in  the  prior  scenarios  and  for  MBL  air  spending  one-twentieth  of  its  time  in  a  cloud  are  shown  in  Fig. 
35.  Solutions  are  shown  for  cases  when  the  H20,  source  SH  are  one-half,  equal  to,  and  double  the  S02 
source  strength.  When  the  H202  source  is  larger  than  the  S02  source  strength,  the  S02  concentrations 
will  be  small. 


TIME  (days) 


Fig.  35  —  Concentrations  of  S02,  H202,  and  aerosol  mass  when  the 
source  rate  of  H202  is  twice  as  large  as  that  for  S02.  If  the  source 
rate  of  H202  is  half  that  of  S02,  then  the  curves  for  S02  and  H202 
are  reversed.  Drshed  lines  are  for  the  case  when  the  source  rate  of 
H202  is  the  same  as  that  of  S02. 
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In  a  more  sophisticated  chemical  model  of  remote  stratiform  clouds,  Chameides  (1984)  examines  S02 
conversion  starting  with  fixed  initial  values  of  S02  and  with  H202  concentrations  generated  by  the 
photochemical  part  of  the  model.  When  the  initial  S02  concentration  was  less  than  about  1  ppb,  sufficient 
H202  was  generated  in  the  standard  stratiform  cloud  model  to  convert  all  SO^  to  sulfate.  Thus  it  would 
appear  that  HiC^  is  more  likely  to  predominate  over  S02  in  the  remote  marine  environment,  resulting 
in  very  low  S02  concentrations.  This  is  in  contrast  to  continental  measurements  made  at  two  mountaintop 
sites  in  the  Eastern  U.S.  (Meagher  et  al.  1990)  where  the  SO^  concentration  was  generally  found  to 
exceed  the  H202. 

It  is  clear  from  the  above  three  scenarios  that  if  either  NH3  or  H202  is  present  in  concentrations 
which  exceed  that  of  S02,  and  if  MBL  clouds  are  present,  then  the  concentration  of  SC^  should  be  very 
low  and  loss  by  gas-phase  and  surface  deposition  will  be  low  compared  to  liquid-phase  conversion.  If 
liquid-phase  oxidation  is  primarily  by  03,  then  S02  concentrations  would  be  higher  than  observed  and 
loss  by  gas-phase  oxidation  and  surface  deposition  would  be  required  to  account  for  the  observed  SO^ 
concentrations.  If  liquid-phase  conversion  is  important,  as  seems  probable,  MBL  cloud  cover  should  be 
an  important  parameter  in  determining  ambient  S02  concentrations. 
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